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ABSTRACT 


EXPERIMENTAL  INVESTIGATION  OF 
HYDRODYNAMIC  SELF-ACTING  GAS  BEARINGS 
AT  HIGH  KNUDSEN  NUMBERS 

This  investigation  gives  experimental  confirmation  of 
the  "slip  flow"  theory,  developed  by  Burgdorfer,  for 
modelling  hydrodynamie  gas  bearing  with  clearances  below 
0.25  pm.  An  inter ferometr ic  technique  using  two  CW  lasers 
is  used  to  measure  the  small  clearances  with  an  0.025  pi 
accuracy.  The  pitch  and  roll  angles  of  the  slider  bearing 
and  the  surface  contour  of  the  slider  are  also  measured 
inter ferometr ically. 

The  effects  of  molecular  rarefaction  are  studied  by 
operating  the  bearing  in  different  gas  media  with 
different  mean  free  paths.  Air  is  used  to  study  Knudsen 
number  below  0.5  while  helium  is  used  to  study  Knudsen 
number  above  0.5. 

Experimentally  measured  trailing  edge  clearances  and 
pitch  angles  are  compared  with  theoretical  predictions 
using  the  "modified"  Reynolds  equation  with  velocity  slip 
boundary  conditions.  The  modified  Reynolds  equation  is 
solved  numerically  using  two  different  numerical 
algorithms.  Excellent  agreement  between  experiment  and 
theory  is  found  for  clearances  as  high  as  1.60  urn  to  as 
low  as  0.075  pm  with  corresponding  ambient  Knudsen  number 
of  0. 04  and  2.51,  respectively. 
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1.  INTRODUCTION 

This  investigation  is  an  experimental  study  of 
hydrodynamic,  self-acting,  gas-lubricated  bearings 
operating  under  "si ip- flow"  conditions.  Conventional 
lubrication  theory  developed  for  gas  bearings  cannot  be 
applied  to  the  type  of  bearing  under  study  because  the 
fluid  inside  the  gas  film  cannot  be  assumed  to  be  a 
continuum,  that  is,  the  molecular  mean  free  path  of  the 
gas  is  not  negligible  compared  with  the  bearing  clearance. 
Under  this  condition,  slippage  between  the  gas  and  the 
bearing  walls  will  occur  and  must  be  considered  in  the 
theoretical  model.  Burgdorfer  [5]  has  developed  a 
slip-flow  model  in  which  first  order  approximation  of  the 
velocity  slip  at  the  boundaries  is  considered  and  a 
so-called  "modified"  Reynolds  equation  is  derived.  This 
equation  is  now  used  as  the  governing  equation  for 
hydrodynamic,  gas-lubricated  bearings  operating  under 
"slip-flow"  conditions.  Burgdorfer  further  asserted  that 
as  long  as  the  ratio  of  the  molecular  mean  free  path  of 
the  gas  and  the  film  thickness  (known  as  the  Knudsen 
number)  is  between  0  and  1,  the  modified  Reynolds  equation 
will  be  sufficient  to  describe  the  flow  phenomenon.  Flows 
with  Knudsen  numbers  between  0.01  and  0.1  are  known  as 
slip  flows,  while  flows  with  Knudsen  numbers  between  0.1 
and  1.0  are  known  as  transition  flows.  Presently  there 
exist  some  experimental  data  to  confirm  the  applicability 
of  the  slip-flow  approximation  to  the  slip  flow  regime. 
However,  there  is  very  little  experimental  data  presently 
to  verify  the  applicability  of  the  slip-flow  approximation 
to  the  transition  flow  regime. 


Experimental  investigations  by  Tseng  [35]  and  Sereny 
[33]  attempted  to  study  the  equation  over  the  full  range 
of  Knudsen  numbers,  that  is,  in  both  slip  flow  and 
transition  flow  regimes.  However,  due  to  the  experimental 
procedure  used  in  their  experiments,  the  effects  of 
molecular  slip  could  not  be  studied  for  Knudsen  numbers 
greater  than  0.1.  In  both  studies,  even  though  the 
ambient  Knudsen  number  was  increased  beyond  0.1  by 
reducing  the  ambient  operating  pressure,  the  additional 
molecular  slippage  effects  expected  as  the  result  of  the 
increase  in  mean  free  path  were  offset  by  the  effects  of 
high  bearing  number  or  compressibility  number,  a  ,  which 
was  also  increased  by  the  lowering  the  ambient  pressure 
since  the  bearing  number  is  inversely  proportional  to  the 
ambient  pressure,  i.e.  A  *  6vUl/ (paht2) .  Thus, 

increasing  the  Knudsen  number  by  decreasing  the  operating 
pressure  does  not  necessarily  increase  the  molecular 
rarefaction  effects. 

Recently,  an  experimental  investigation  by  Kogure, 
Kaneko,  and  Hitsuya  [21]  has  shown  that  the  modified 
Reynolds  equation  is  valid  for  Knudsen  number  up  to  0.3. 
In  their  work,  the  Knudsen  number  was  increased  by 
decreasing  the  bearing  clearance  through  reducing  the 
bearing  width.  In  contrast  to  the  work  of  Tseng  [35]  and 
Sereny  1 33  D  *  the  increased  molecular  slip  effect  in 
Kogure 's  work  was  not  overshadowed  by  the  high  bearing 
number  effect.  However,  this  recent  study  by  Kogure  was 
limited  to  the  low  end  of  the  transition  flow  regime. 

In  the  present  experimental  investigation,  in  order 
to  minimize  the  high  bearing  number  effects  and  maximize 


the  molecular  rarefaction  effects,  the  bearing  is  tested 
in  different  gas  media  of  different  ambient  mean  free 
paths.  More  specifically,  air  is  used  to  study  Knudsen 
numbers  below  0.5  while  helium  is  used  to  study  Knudsen 
numbers  above  0.5.  Furthermore,  slider  bearings  of 
various  width  are  used  to  change  the  bearing  clearance 
which  in  turn  will  change  the  Knudsen  number.  This 
investigation  represents  the  first  work  that  has 
successfully  studied  the  hydrodynamic,  self-acting, 
gas-lubricated  bearing  operating  over  a  wide  range  of 
Knudsen  numbers  (that  is,  in  both  slip  flow  and  transition 
flow  regimes)  without  the  presence  of  high  bearing  number 
effects . 

The  experimental  results  of  this  investigation 
confirm  the  validity  of  the  modified  Reynolds  equation  for 
the  entire  reported  range  of  Knudsen  numbers,  between  0 
and  1.  Moreover,  the  results  also  showed  that  the 
modified  Reynolds  equation  is  also  capable  of  describing 
the  flow  phenomenon  for  ambient  Knudsen  number  far  beyond 
1.0. 

In  addition  to  the  experimental  work,  the 
investigation  also  examined  the  accuracy  of  the  numerical 
solution  of  the  modified  Reynolds  equation.  Two 
finite-difference  algorithms  were  used  to  solve  the 
equation.  Numerical  results  -  the  predicted  load  and 
pitch  angle  -  from  the  two  algorithms  agree  with  each 
other  to  within  2%. 


2.  STATE  OF  THE  ART 

The  study  of  ultra  thin  gas  film  bearing  has  become 
of  great  interest  in  recent  years.  One  of  the  prime 
motivations  has  come  from  the  computer  magnetic  disk 
recording  industry  where  higher  recording  density  and 
signal  resolution  can  be  achieved  if  the  read/write 
element  which  is  attached  to  a  taper-flat  slider  can  be 
maintained  at  a  closer  spacing  over  the  magnetic  disk 
rotating  at  a  high  speed.  It  has  been  established  by 
Hoagland  [19]  and  Pear  [24]  that  the  recording  density  can 
be  greatly  increased  if  the  slider/disk  separation  is 
reduced.  In  actual  operation,  due  to  the  relative  motion 
of  the  disk  and  the  taper-flat  slider,  a  hydrodynamic 
self-acting  gas  lubricated  bearing  is  established  between 
the  two  surfaces.  The  total  pressure  force  generated 
within  the  gas  film  balances  exactly  the  externally 
applied  force  on  the  slider.  The  physical  situation 
described  is  illustrated  schematically  in  Figure  11-2. 
Thus,  to  design  a  recording  system  where  a  proper 
slider/disk  separation  can  be  maintained,  an  accurate 
model  of  the  bearing  system  is  needed. 

Since  the  first  conception  and  application  of  the 
disk  recording  head  element  in  1956,  the  classical 

Reynolds  equation  has  been  found  to  predict  well  the 
performance  of  the  air  bearing  operating  with  film 

thicknesses  on  the  order  of  8  to  10  ym  [4].  The  sliders 
studied  by  [4]  were  either  plane  flat  sliders  or  convex 

sliders  with  very  large  width  to  length  ratio  or 

slenderness  ratio  (about  1.3) •  As  the  slider  underwent 
design  changes  in  order  to  increase  the  recording  density, 


the  result  was  a  reduction  in  the  bearing  clearance  from  8 
Vin  to  1  ym  [35].  The  single  pad  slider  geometry  also 
became  more  sophisticated  with  a  taper-flat  surface  and  a 
slenderness  ratio  of  2.0.  To  further  increase  the 
recording  density,  the  bearing  clearance  was  further 
reduced  from  1  yra  to  the  present  0.5  ym.  The  slider 
geometry  was  also  changed  to  the  present  double  pad 
arrangement  with  taper-flat  surfaces  and  a  slenderness 
ratio  of  0.1.  This  type  of  head  is  known  as  the 
Winchester  head  and  is  shown  schematically  in  Figure  11-3. 

As  the  bearing  clearance  is  reduced  to  the  present 
0.5  ym  value,  the  classical  Reynolds  equation  used  to 
describe  the  fluid  flow  phenomenon  in  the  bearing  film  is 
no  longer  adequate.  It  has  been  shown  by  Burgdorfer  [5] 
that  as  the  slider/disk  separation  is  substantially 
reduced,  the  classical  Reynolds  equation  must  be  modified 
to  take  into  account  molecular  slip  at  the  gas/solid 
interfaces.  The  main  effect  of  this  molecular  rarefaction 
is  to  reduce  the  load  carrying  capacity  of  the  slider 
bearing.  Figure  11-4  shows  the  extent  of  this  molecular 
slip  effect  as  the  mean  free  path  of  the  gas  medium  is 
increased  at  constant  ambient  pressure.  As  the  figure 
indicates,  for  fixed  trailing  edge  clearance,  less 
external  load  is  required  to  maintain  the  bearing  at  that 
clearance  when  slip  is  present.  In  other  words,  if  the 
load  on  the  slider  is  fixed,  the  slider  will  operate  at  a 
lower  clearance  when  molecular  slip  is  considered.  Thus, 
if  molecular  slip  is  not  properly  accounted  for  in 
designing  the  bearing  system,  substantial  error  can  result 
in  predicting  the  bearing  clearance  for  constant  load. 


The  modified  Reynolds  equation  is  theoretically  valid  for 
Knudsen  number  up  to  1.  However,  there  exist  very  little 
experimental  data  presently  to  substantiate  the  slip-flow 
approximation  of  the  Reynolds  equation  developed  by 
Burgdorfer  [5]  over  the  entire  range  of  Knudsen  numbers 
( from  0  to  1.0). 

Hsing  and  Malanoski  [20]  were  the  first  investigators 
to  attempt  to  confirm  the  validity  of  the  "slip- flow" 
theory  as  applied  to  a  spiral  grooved  thrust  bearing. 
They  compared  their  theoretically  predicted  pressures  to 
those  measured  by  Sternak  [20]  and  found  significant 
improvement  in  the  theoretically  results  when  slip  was 
considered.  Excellent  agreement  was  found  for  certain 
range  of  operating  conditions  -  compressibility  number 
less  than  100  and  ambient  Knudsen  number  less  than  0.1. 
However,  for  conditions  above  this  range,  the  slip-flow 
theory  begins  to  fail  quite  rapidly.  Consequently,  Hsing 
and  Malanoski  seriously  questioned  the  validity  of  the 
slip-flow  approximation  when  the  ambient  Knudsen  number 
becomes  greater  than  0.1,  that  is  in  the  transition  flow 
regime . 

The  first  experimental  work  that  attempted  to  study 
the  slip-flow  approximation  as  applied  to  the  slider 
bearing  system  was  carried  out  by  R.C.  Tseng  [35].  In  his 
work,  Tseng  was  interested  primarily  in  the  molecular  slip 
effects  in  an  "infinitely"  wide  bearing.  Thus,  in  his 
experiments,  a  taper-flat  slider  with  a  very  large  width 
to  length  ratio  (about  2.0)  was  used.  To  study  the 
rarefaction  effects,  the  molecular  mean  free  path  of  the 
air  was  increased  by  decreasing  the  ambient  pressure.  In 


this  manner,  Tseng  was  able  to  run  the  wide  slider  bearing 
with  ambient  Knudsen  numbers  up  to  2.4.  It  should  be 
noted  that  this  is  far  beyond  the  reported  limits  of  the 
slip  flow  approximation.  However,  even  though  such  high 
ambient  Knudsen  numbers  were  achieved,  the  effect  of  slip 
was  negligible  because  the  bearing  performance  was 
dominated  by  high  compressibility  (bearing)  number  effect. 
This  is  particularly  evident  when  one  examines  Tseng’s 
figure  since  his  continuum  theory  curve  and  the  slip-flow 
theory  curve  coincide  at  the  high  Knudsen  number  range. 
Nevertheless,  Tseng  was  able  to  confirm  that  the  slip  flow 
approximation  can  predict  the  trailing  edge  clearance  to 
within  10%  of  the  actual  value  for  Knudsen  numbers  up  to 
0.1  (slip  flow  regime).  Tseng's  study  also  found  that  a 
proper  selection  of  the  surface  accommodation  coefficient 
for  air/disk  and  air/slider  interfaces  can  greatly  improve 
the  theoretical  predictions.  Surface  accommodation 
coefficient,  a  ,  is  defined  as  the  fraction  of  gas 
molecules  which  are  adsorbed  when  striking  the  wall  and 
then  reemitted  with  a  random  distribution  in  angle  and 
velocity.  Tseng  found  that  if  the  surface  accommodation 
coefficient  is  taken  to  be  0.89  instead  of  1.0,  a  closer 
agreement  between  theory  and  experiment  can  be  achieved. 
It  should  be  noted  that  the  value  0.89  used  for  the 
surface  accommodation  coefficient  was  the  experimental 
value  measured  by  Schaaf  £303  for  an  air/glass  interface 
for  flows  with  low  shear  rates.  In  the  bearing  film,  the 
shear  rate  is  several  orders  of  magnitude  larger  than 
those  encountered  by  Schaaf  in  his  studies.  Whether  the 
same  value  can  be  used  is  debatable.  In  the  present  study, 
the  experiments  have  been  used  to  help  in  the 
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determination  of  the  best  value  for  the  surface 
accommodation  coefficient  for  air/glass  and  helium/glass 
interfaces . 

Only  recently  have  experimental  data  become  available 

for  narrow  slider  bearings  with  width  to  length  ratio  of 

0.1.  The  data  were  obtained  by  Sereny  [333-  With  a 

narrower  bearing,  the  effect  of  slip  is  expected  to  be 

more  pronounced.  Similar  to  Tseng's  work  [35],  Sereny 

studied  the  molecular  rarefaction  effect  by  increasing  the 

mean  free  path  of  air  via  decreasing  the  ambient  pressure. 

In  his  study,  the  effect  of  surface  curvature  (crown)  due 

to  lapping  was  considered  because  it  was  found  that  the 

crown  could  have  considerable  effect  on  the  bearing 

performance.  Figure  11-5  is  a  typical  curve  showing  the 

effects  of  crown  on  the  load  and  pitch  angle  calculations. 

As  one  can  see,  a  convex  crown  (positive  hc)  can 

significantly  increase  the  load  capacity  while  a  concave 

C  crown  (negative  hc)  decreases  the  load  capacity.  Based  on 

his  experimental  results,  Sereny  concluded  that  the 

slip-flow  approximation  completely  fails  to  predict  the 

rarefaction  effects  for  ambient  Knudsen  number  greater 
t 

than  0.1.  It  should  be  noted  that  theoretical 

load/spacing  curves  plotted  did  not  shift  as  the  ambient 
operating  pressure  is  reduced  at  constant  disk  velocity. 
^  This  is  not  surprising  since  as  the  Knudsen  number  is 

increased  by  decreasing  the  ambient  pressure,  the 

operating  bearing  number  is  also  increased  by  the  same 
t  factor.  Consequently,  the  additional  molecular  slip 

effect  as  the  result  of  the  increase  in  mean  free  path  is 
offset  by  the  increased  high  bearing  number  effect. 


t 


Figure  11-6  shows  this  offsetting  effect  which  is 
illustrated  by  the  fact  that  the  slip  flow  curves  for 
subambient  pressures  do  not  shift  away  from  ambient 
pressure  curve.  Comparing  figure  11-4  with  figure  11-6, 
it  is  clear  that  by  changing  the  mean  free  path  while 
keeping  the  ambient  pressure  constant,  more  slippage 
effect  could  be  observed. 

Recently  Kogure,  Kaneko ,  and  Mitsuya  [21]  reported 
their  experimental  findings  on  narrow  slider  bearings 
operating  at  "high"  Knudsen  numbers.  They  found  excellent 
agreement  between  experiment  and  theory  for  Knudsen 
numbers  up  to  0.35.  Unlike  previous  experimental 
investigations  [33,  35],  Kogure  studied  the  increased 
molecular  slip  effect  by  decreasing  the  trailing  edge 
clearance  via  reduction  of  the  bearing  width.  Figure  11-7 
is  a  typical  curve  showing  the  sensitivity  of  the  load 
capacity  to  bearing  width.  In  this  manner,  high  bearing 
numbers  encountered  by  Tseng  [35]  and  Sereny  [33]  could  be 
avoided.  Kogure' s  study  was  limited  to  Knudsen  numbers 
below  0.35  because  the  mean  free  path  of  air  (0.069  inn) 
was  not  changed  and  the  lowest  trailing  edge  clearance 
achieved  was  about  0.2  ym.  However,  for  Knudsen  number 
range  greater  than  0.35,  there  is  still  no  experimental 
data  available  to  either  prove  or  disprove  the  slip-flow 
approximation . 

2.1  Measurement  Techniques 

Capacitance  and  optical  techniques  have  been  used 
successfully  to  measure  the  slider/disk  separation. 
Brunner  and  his  colleagues  [4]  used  two  capacitance  probes 
imbedded  in  the  slider  to  measure  the  bearing  clearances 


at  two  points  on  the  slider.  Knowning  the  surface 

contour,  Brunner  was  able  to  determine  the  trailing  edge 
clearance  as  well  as  the  pitch  angle  of  the  slider 

bearing.  It  is  not  known  whether  the  probes  had  any 

effect  on  the  bearing  performance.  As  the  slider  size  was 
reduced,  smaller  and  unshielded  capacitance  probes  were 
used  [3]  to  measure  the  slider  flying  characteristics. 
However,  with  the  present  slider  design,  neither  the 

shielded  nor  the  unshielded  could  be  used  because  the 
probes  are  much  larger  than  the  slider  itself. 

Optical  interferometric  techniques  have  been 
successfully  used  to  measure  the  gap  between  a  moving 
magnetic  tape  and  a  static  recording  head  [23,  18].  The 
interference  fringe  patterns  can  be  produce  by  either  a 
monochromatic  light  source  or  a  white  light.  Lin  and 
Sullivan  [22]  were  the  first  investigators  to  successfully 
use  white  light  to  produce  color  fringes  (known  as  Newton 
color  rings)  for  measuring  the  trailing  edge  clearance  of 
slider  bearings.  Later  Tseng  [35]  also  chose  to  use  white 
light  over  monochromatic  light  to  measure  the  slider 
trailing  edge  clearance.  Recently  Sereny  [33)  used  two 
lasers  (a  CW  He-Ne  laser  and  a  Variable  Wavelength  Pulsed 
Dye  laser)  to  determine  the  clearance  profile  of  the 
slider  bearing.  Since  the  optical  measuring  techniques 
gives  a  complete  mapping  of  the  slider  clearance,  Sereny 
also  used  this  technique  to  measure  the  surface  curvature 
(crown)  on  the  "flat"  portion  of  the  taper-flat  slider  as 
well  as  the  taper  height. 

There  are  advantages  and  disadvantages  associated 
with  the  two  techniques.  For  capacitance  probes,  the 


advantage  is  that  a  real  disk  can  be  used  in  the 
experiments,  while  the  disadvantage  is  that  the  imbedded 
probe  can  interfere  with  the  bearing  performance  due  to 
added  weight  or  due  to  improper  probe  implantation.  On 
the  other  hand,  the  advantage  of  the  interferometric 
technique  is  that  the  slider  is  not  altered  in  any  way 
which  could  interfere  with  its  performance,  while  the 
disadvantage  is  that  the  real  disk  must  be  replaced  with  a 
glass  disk  in  order  to  allow  the  light  to  pass  through  and 
to  creat  the  necessary  optical  interference.  Slider 
flying  characteristics  can  be  different  if  the  glass  disk 
is  used  because  the  surface  finish  (roughness)  on  an 
actual  disk  is  different  from  the  glass  disk.  However,  it 
is  well  accepted  now  that  the  inter ferometr ic  technique  is 
a  reliable  and  accurate  method  for  measuring  small 
clearances . 

2.2  Surface  Roughness 

As  the  bearing  clearance  is  reduced,  effects  due  to 
surface  roughness  on  the  slider  and  on  the  disk  may  become 
significant  in  addition  to  the  molecular  rarefaction 
effects.  Roughness  studies  by  numerous  researchers 
[7,  14,  27,  34,  36]  have  shown  that  surface  roughness  can 
considerably  influence  the  bearing  performance.  The 
roughness  can  be  transverse,  longitudinal,  or  striated 
relative  to  the  direction  of  the  flow.  Rhow  and  Elrod 
[27]  have  shown  that  there  is  a  general  increase  in  the 
load  carrying  capacity  due  to  striated  roughness  based  on 
mean  film  thickness.  However,  a  theoretical  study  based 
on  Reynolds  roughness  by  Christensen  and  Tonder  [7]  has 
shown  that  effects  due  to  transverse  or  longitudinal 


roughness  are  less  pronounced  in  the  case  of  narrow 
bearings  (width  to  length  ratio  of  0.5)  as  compared  to 
"infinitely"  wide  bearing.  In  the  present  study,  sliders 
with  width  to  length  ratios  of  0.15  or  less  are  studied 
and  surface  roughness  is  not  considered. 


In  the  present  experimental  study,  the  molecular 
rarefaction  effect  predicted  using  the  slip-flow 
approximation  of  the  Reynolds  equation  is  verified 
experimentally.  Bearing  clearance  profiles  and  surface 
contour  are  measured  inter ferometr ically  using  two  CW 
lasers  -  a  Helium-Neon  laser  with  wavelength  of  632.8  nm 
and  a  Helium-Cadmium  laser  with  wavelength  of  441.6  nm. 
Experiments  are  performed  with  ambient  Knud  sen  number 
varying  from  0.04  to  2.5  without  approaching  excessively 
high  bearing  numbers.  Bearing  clearances  as  low  as  0.075 
pm  with  ambient  mean  free  path  of  0.188  pm  have  been 
measured.  Excellent  agreement  has  been  found  between 
experiment  and  slip-flow  theory  over  the  entire  range  of 
ambient  Knudsen  number,  from  0.04  to  2.5. 


3.  THEORY 

The  classical  treatment  of  the  slider  gas  bearing 
problem  is  well  known  and  has  been  thoroughly  discussed  in 
many  lubrication  texts  C  9  >  15,  16].  The  classical 
Reynolds  equation  is  limited ,  however,  to  flows  where  the 
molecular  mean  free  path  of  the  gas  is  negligible  compared 
with  the  dimensions  of  the  flow  passage,  namely  the  slider 
bearing  clearance.  On  the  other  hand,  when  the  film 
thickness  is  much  smaller  than  the  molecular  mean  free 
path  of  the  gas  (i.e.  rarefied  gas  flow),  the  problem  is 
solved  using  kinetic  theory  analysis.  However,  if  the 
flow  is  neither  continuum  nor  molecular,  the  classical 
treatment  is  inadequate  and  kinetic  theory  analysis  is  too 
tedious.  It  has  been  shown  from  rarefied  gas  studies, 
both  theoretically  and  experimentally,  that  the 
Navier-Stokes  equation,  together  with  slip  velocity  and 
temperature  jump  conditions,  is  adequate  for  analyzing 
flows  in  the  so-called  "slip-flow”  or  "transition"  regime. 

The  first  theoretical  analysis  of  the  rarefaction 
effect  in  gas  lubrication  was  done  by  Burgdorfer  [5].  The 
derivation  of  the  modified  Reynolds  equation  is  identical 
to  that  for  the  classical  Reynolds  equation  except  for  the 
boundary  conditions  at  the  walls  (i.e.  the  slider  surface 
and  the  rotating  disk  surface.)  In  addition  to  the 
standard  Couette  flow  boundary  conditions,  a  velocity  slip 
is  also  assumed  at  all  surface  boundaries.  The  magnitude 
of  the  slip  is  approximated  to  the  first  order  by  the 
expression:  (refer  to  coordinate  system  defined  in  Figure 
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f  -  a  numerical  constant 
a  -  surface  accommodation  coefficient 
X  -  mean  free  path  of  gas 


In  equation  (3.1),  f  is  usually  taken  to  be  unity  while 
a  can  be  less  than  or  equal  to  unity,  depending  on  the  gas 
and  boundary  surface  finish.  No  temperature  jump 
condition  is  needed  since  theoretical  estimates  of 
temperature  difference  between  the  bearing  surface  and  the 
gas  film  show  this  to  be  negligible. 

For  the  sake  of  simplicity,  the  modified  Reynolds 
equation  is  usually  written  with  o  and  f  equal  to  1.  But 
if  either  parameter  is  not  unity,  the  equation  can  be 
easily  modified  by  introducing  the  factor 
f (2  - o  )/ a  wherever  appears.  The  following  steady-state 
modified  Reynolds  equation  was  derived  by  (53: 


&03H<t  +  £K!?('  +tM  -  6>,u^(ph)  u.2> 

Nondimensionalizing  the  variables  as  defined  in  the 
Nomenclature,  the  resulting  dimensionless  form  of  equation 
(3.2)  is  written  as  follows: 

Cr)’MPH3!y(’  ♦  Ttr>]  +  *  Tff']-  a(t)2 ^r1 

(3.3) 


where  K„is  the  ratio  of  the  ambient 


■•an  free  path  to  the  minimum 
baaring  claaranca  (known  as  the 
ambient  Knud sen  number) 

Equation  (3.3)  is  written  in  slightly  different  form 
than  the  equation  derived  by  [5]  in  that  the  ratio  of  the 
slider  width  to  the  slider  length  is  present.  Instead  of 
the  normal  procedure  of  normalizing  the  y-dimension  with 
respect  to  the  length,  normalization  is  carried  out  with 
respect  to  the  width.  For  sake  of  completeness,  the 
derivation  steps  leading  to  equation  (3. 3)  are  outlined  in 
Appendix  A. 

If  o  takes  on  values  less  than  unity,  equation  (3*3) 

is  slightly  modified  with  the  parameter  being  replaced 

by  (2  -  o)K  / a  .  For  air/glass  interface,  it  was  found 
"  00 

by  Schaaf  t30j  that  a  realistic  value  of  o  is  0.89. 

Moreover,  Tseng  [35]  and  Sereny  [33]  both  have  shown  that 

with  o  equal  to  0.89,  a  better  agreement  between 
experiment  and  slip-flow  theory  can  be  achieved. 

Equation  (3.3)  is  a  two-dimensional,  non-linear, 
second-order  differential  equation  in  P  and  cannot  be 
solved  analytically  for  arbitrary  slider  geometry.  But 

several  observations  can  be  made  by  examining  the 
equation.  One  is  that  if  the  bearing  number,  A  ,  is 

sufficiently  large  and  the  width  to  length  ratio  is  of 
order  1,  the  right  hand  side  of  the  equation  will 
dominate.  The  flow  in  the  gas  film  is  therefore  Couette 
dominant  and  molecular  slip  effects  are  negligible. 
Another  observation  is  that  if  the  slider  is  very  narrow 
(width  to  length  ratio  ^  0.1),  as  is  the  case  for  the 
present  bearing  geometry,  the  flow  out  the  side  edges  will 


be  significant.  Moreover,  the  criterion  for  "large" 
bearing  number  is  now  dictated  by  the  new  dimensionless 
parameter,  a*  (  =  A  <w/l)2)  ,  the  "modified"  bearing 

number.  Thus  for  narrow  bearings,  the  flow  inside  the  gas 
film  will  not  be  Couette  dominant  until  the  modified 
bearing  number,  a*,  is  sufficiently  large.  Therefore, 
until  this  limit  is  reached,  the  molecular  slip  effect  in 
narrow  bearings  cannot  be  neglected  even  though  the 
conventional  bearing  number  is  high. 

The  ability  to  obtain  analytical  solutions  to  any 
problem  is  important  because  it  allows  parametric  analyses 
which  can  lead  to  better  understanding  of  the  problem. 
Simplification  of  the  governing  equation  is  one  of  the 
ways  that  can  lead  to  analytical  solutions.  In  gas 
bearings,  the  most  common  way  to  simplify  equation  (3.3) 
is  to  assume  the  bearing  is  infinitely  wide.  With  this 
assumption,  the  equation  becomes  one-dimensional  and  can 
be  easily  solved  for  certain  slider  geometries. 
Burgdorfer  was  able  to  obtain  closed  form  solutions  for 
Rayleigh  step  bearings  and  plane  slider  bearings  [53. 
Perturbation  techniques  and  asymptotic  analysis  have  been 
used  by  various  investigators  to  gain  better  insight  into 
the  problem.  Asymptotic  solution  for  infinitely  wide 
slider  bearings  at  high  bearing  numbers  without  slip  have 
been  obtained  by  numerous  investigators 

t  * 

[10,  11,  13,  17,  313-  Recently,  Sereny  obtained  the 

one-dimensional  asymptotic  solution  with  slip  at  high 
bearing  numbers  [333.  Asymptotic  solutions  for  finite 
width  slider  bearings  at  high  bearing  numbers  without  slip 
have  been  obtained  by  Elrod  and  McCabe  [133  and  Schmitt 
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and  DiPrima  [32].  However,  an  asymptotic  solution  for 
finite-width  slider  bearings  at  high  bearing  numbers  with 
slip  has  not  been  developed  thus  far. 

In  order  to  design  and  develop  new  bearing  systems, 
the  full  modified  Reynolds  equation,  eq.  (3.3) »  must  be 
considered  and  the  solution  can  be  obtained  only  by 
numerical  techniques. 
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4.  NUMERICAL  SOLUTION 

The  analysis  of  compressible  lubrication  films  is 
complicated  due  to  the  non-linear  nature  of  the  governing 
equation,  the  well  known  Reynolds  equation.  The  equation 
is  further  complicated  by  the  introduction  of  molecular 
slip  as  the  bearing  clearance  become  comparable  to  the 
molecular  mean  free  path  of  the  gas  medium.  Analytical 
solution  of  the  full  modified  Reynolds  equation  is 
impossible  to  obtain  except  for  certain  limiting  cases 
discussed  in  the  previous  chapter.  For  real  problems,  the 
full  equation  must  be  considered  and  solved  numerically. 


Presently,  there  are  two  schools  of  thought  on  the 
type  of  numerical  solution  technique  that  should  be  used. 
The  majority  of  the  numerical  solution  are  still  based  on 
the  finite  difference  approximation  of  the  governing 
equation.  But  in  recent  years,  the  finite  element  method, 
whose  application  had  been  limited  mainly  to  structural 
analysis,  has  been  applied  to  lubrication  problems  (as 
well  as  other  fluids  and  heat  transfer  problems.) 

Several  papers  have  been  published  recently  which 
focused  on  the  application  of  the  finite  element  technique 
in  solving  the  Reynolds  equation  [2,  25,  26).  The  finite 
element  method  has  one  distinct  advantage  over  the  more 
commonly  used  finite  difference  method;  as  the  former  is 
able  to  handle  irregular  shape  boundaries  without  major 
modifications  to  the  computer  code.  This  versatility  is 
attributed  to  the  fact  that  irregular  shapes  can  always  be 
approximated  by  patch  work  of  polygonal  elements. 
However,  both  Booker  and  Huebner  [21  and  Reddi  and  Chu 
[26]  are  careful  to  point  out  that,  in  terms  of  computer 


tine,  the  finite  element  method  may  not  be  better 
(sometimes  even  worse)  than  the  finite  difference  method. 
Reddi  and  Chu  [26]  also  noted  that  it  is  very  difficult  to 
compare  the  two  techniques  on  the  same  level  since  the 
solution  techniques  for  solving  the  resulting  matrices  are 
different . 

The  most  often  used  numerical  technique  to  solve  the 
full  compressible  Reynolds  equation  is  still  the  finite 
difference  method.  Castelli  and  Pirvics  [6]  summarized  in 
their  paper  the  state  of  the  art  of  numerical  solution 
methods  as  it  existed  in  1969.  Since  that  time  only  a  few 
papers  [8,  37.  12]  have  focused  on  the  numerical  solution 
techniques  as  applied  to  the  Reynolds  equation.  The 
finite  difference  method  is  divided  into  two  distinct 
parts  -  discretization  of  the  Reynolds  equation  and  the 
solution  of  the  resulting  set  of  algebraic  equations. 

All  numerical  solution  methods  require  the 
representation  of  the  dependent  variable  at  a  finite 
number  of  points  defined  by  the  intersection  of  a  grid 
mesh.  The  distribution  of  grid  lines  should  be  such  that 
they  fall  on  all  boundaries  and  clearance  discontinuities. 
Furthermore,  if  there  are  sharp  pressure  gradients 
anywhere  in  the  bearing,  special  grid  distribution  is 
required.  For  the  taper-flat  slider  bearing,  sharp 
pressure  gradients  exist  at  the  break  and  at  the  trailing 
edge  and  denser  grid  mesh  is  need  at  these  two  regions  in 
order  to  assure  accurate  solution  and  numerical  stability. 
These  sharp  gradients  are  caused  by  the  presence  of  small 
boundary  layers  and  these  boundary  layers  have  been 
estimated  to  be  on  the  order  of  1/A  [10,  11,  13,  31,  32]. 


For  finite  width  bearings,  if  the  bearing  number  is 
sufficiently  large,  a  side  boundary  layer  of  the  order 
1//A  C 1 3 »  32]  will  be  present  at  the  edges  causing  a 
large  pressure  gradient  in  the  region  and  a  finer  grid 
mesh  will  be  required  at  the  side  edges  of  the  slider  for 
numerical  stability. 

Until  recently,  all  finite  difference  algorithms  that 
have  been  developed  [6,  37  3  require  the  use  of  variable 
grid  mesh  when  the  bearing  number  is  large.  However,  a 
new  algorithm  has  been  developed  by  H.G.  Elrod  [12]  which 
does  not  require  the  use  of  variable  grid  spacing  in  order 
to  insure  numerical  stability.  Because  constant  grid  mesh 
can  be  used,  this  new  algorithm  is  much  faster  than  any 
previous  algorithm  in  terms  of  computer  time.  The 
original  algorithm  developed  by  Elrod  has  not  been 
published.  In  this  investigation  a  similar  algorithm  was 
formulated  using  the  exponential  approximation  method 
developed  by  Elrod  [12].  The  formulation  of  this  "new" 
algorithm  is  outlined  in  Appendix  B. 

In  this  study,  two  different  discretization  methods 
are  used.  One  is  the  integral  discretization  technique 
developed  by  Castelli  and  Pirvics  [6].  The  details  of 
this  technique  can  be  found  in  Reference  [6],  Basically 
for  integral  discretization,  a  finite  difference 
approximation  is  applied  to  equation  (3*3)  after  the 
equation  has  been  integrated  once  over  each  grid  cell 
using  Gauss'  theorem  (for  transforming  surface  integral  to 
line  integral).  It  was  shown  by  Sereny  [333  that  integral 
discretization  is  superior  to  the  differential 
discretization  in  that  it  is  less  dependent  on  the  size  of 


the  grid  mesh  and  provides  the  necessary  conservation  of 
mass.  The  other  discretization  method  used  is  the 
exponential  approximation  method  developed  by  Elrod  [12]. 

Once  the  Reynolds  equation  is  discretized,  one 
algebraic  equation  is  written  for  each  grid  point 
resulting  in  a  system  of  algebraic  equations.  The  system 
can  then  be  solved  by  various  numerical  techniques.  Two 
implicit  schemes  are  used  -  Column  Method  and  Alternating 
Directional  Implicit  method  (ADI).  The  so-called  Column 
Method  is  described  in  detail  in  reference  [6]. 

In  the  beginning  of  this  investigation,  experimental 
results  are  compared  with  theoretical  predictions  using 
the  integral  discretization  techniques.  But  later  the 
exponential  approximation  method  is  used  for  all  the 
predictions.  Both  Column  Method  and  ADI  were  used  with 
each  discretization  method  and  their  results  checked  with 
each  other.  In  addition,  the  efficiency  in  terms  of 
computational  time  of  each  algorithm  are  also  compared. 


5.  DESCRIPTION  OF  THE  EXPERIMENTAL  APPARATUS 

The  measurement  of  film  thicknesses  on  the  order  of 
0.25  yn  or  less  requires  a  clean  and  vibration  free 
environment.  In  addition,  close  tolerances  on  the 

measured  quantities  -  bearing  clearance  profile  and  load  - 
must  be  maintained.  An  experimental  apparatus  meeting 
such  requirements  is  used  in  the  present  investigation. 
The  apparatus  was  originally  designed  and  built  by  Sereny 
C  33  3 .  Two  major  changes  have  been  made  to  the 

experimental  set-up.  First  is  that  the  Variable 

Wavelength  Pulsed  Dye  Laser  has  been  replaced  with  a 
Continuous  Wavelength  (CW)  Helium-Cadmium  Laser.  It  was 
found  that  the  Pulsed  Dye  laser  was  too  unreliable  and  the 
lasing  wavelength  of  the  dye  solution  used  in  the  laser 
could  not  be  determined  accuratetly.  The  accuracy  in 
determining  the  bearing  clearance  profile  is  highly 
dependent  on  knowing  the  wavelength  of  the  two  lasers 
used.  Second  is  that  the  test  chamber  has  been 

instrumented  with  a  gas  analyzer  to  measure  the 

concentration  of  different  gas  media  which  may  be 
introduced  into  the  chamber.  Only  gases  with  sufficiently 
different  thermal  conductivity  than  air  can  be  used. 

The  experimental  apparatus  consists  of  two  major 
components  -  a  sealable  but  controllable  test  chamber  and 
an  optical  bench.  Figure  11-8  is  a  schematic  diagram  of 
the  experimental  set-up. 

5.1  Test  Chamber 

The  test  chamber  is  divided  into  upper  and  lower 
sections  by  a  dividing  plate.  The  variable  speed  D.C. 
motor  and  the  belt  drive  mechanism  are  housed  in  the  lower 


chamber  while  the  actual  bearing  test  environment  is 
located  in  the  upper  chamber.  This  is  done  to  prevent 
contaminants  from  the  motor  and  belt  drive  from  entering 
the  test  area.  The  entire  test  chamber  can  be  sealed  so 
that  the  inside  environment  can  be  controlled  such  that 
the  mean  free  path  can  be  changed.  There  are  two  ways  to 
change  the  mean  free  path  inside  the  test  chamber.  The 
first  method  is  to  reduce  the  chamber  pressure  since  the 
mean  free  path  of  a  gas  is  inversely  proportional  to  the 
pressure.  Alternatively,  the  gas  inside  the  chamber  can 
be  changed  to  another  one  that  has  a  different  mean  free 
path.  For  example,  helium  has  a  mean  free  path  that  is  2.5 
times  greater  than  air  while  the  viscosity  of  the  two 
gases  differ  only  by  ten  percent.  In  the  present 
investigation,  even  though  the  apparatus  was  originally 
designed  to  operate  at  subambient  pressures,  the  second 
method  is  used  because  the  molecular  slip  effect  is  more 
pronounced  if  the  mean  free  path  can  be  increased  without 
changing  the  ambient  pressure.  Comparison  of  Figures  11-4 
and  11-6  clearly  shows  the  superiority  of  the  second 
method.  It  should  be  noted  that  even  though  in  Figure 
11-4  the  curves  are  for  air,  it  is  still  applicable  for 
helium  because  the  viscosity  of  the  two  gases  do  not 
differ  significantly. 

The  upper  and  lower  chambers  are  connected  only 
through  a  valve.  When  the  test  chamber  is  sealed  for 
helium  experiments,  this  valve  is  opened  to  allow  helium, 
which  is  forced  fed  into  the  upper  chamber,  to  flow  into 
the  lower  chamber  before  exiting  the  test  chamber.  In 
this  manner,  there  is  always  a  flow  from  upper  to  lower 


chamber  which  helps  to  prevent  dust  in  the  lower  chamber 
from  migrating  up  to  the  upper  chamber. 

In  the  upper  chamber,  a  4.45  cm  thick  optical-flat 
disk  with  a  25.4  cm  O.D.  is  mounted  at  its  inner  diameter 
(4.6  cm  I.D.)  to  a  precision  ball  bearing  spindle.  The 
disk  has  a  surface  flatness  of  1/20  of  the  He-Ne  laser 
wavelength  of  light  used  and  is  dynamically  balanced  to 
within  2.5  pm  peak  to  peak. 

The  slider  bearing  used  in  this  study  is  of  the  IBM 
"Winchester"  design.  Basically,  the  slider  has  a  two-rail 
taper-flat  configuration  and  is  made  of  ferrite  material. 
A  sketch  of  the  Winchester  slider  is  shown  in  Figure  11-3. 
The  nominal  length  associated  with  the  slider  of  this  type 
is  5.6  mm  with  a  projected  taper  length  of  1.27  mm.  The 
nominal  taper  angle  is  about  8.7  mill i-radians  (half  a 
degree).  In  this  study,  the  nominal  width  of  the  rail 
varies  between  0.38  ram  and  0.89  mm. 

The  slider  head  is  aligned  and  glued  onto  a  gimbal 
spring  suspension  system  which  is  capable  of  applying  a  10 
gram  load  on  the  slider.  The  load  is  applied  through  a 
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load  point  called  the  pivot  point.  The  spring  is  designed 
such  that  very  little  moment  is  created  about  the  pivot 
point  when  the  slider  pivots.  Figure  11-9  is  a  photograph 
showing  the  gimbal  spring  with  the  slider  mounted.  The 
slider-spring  assembly  is  clamped  to  one  end  of  the 
cantilever  beam  while  the  other  end  of  the  beam  is  fixed 
onto  a  precision  stage  with  five  degrees  of  freedom  - 
three  rotation  and  two  translation.  With  this  mounting 
scheme,  the  yaw  and  roll  angles  of  the  slider  can  be 


adjusted  or  biased 


The  external  load  on  the  slider  can  be  varied  through 
the  use  of  a  variable  loading  mechanism.  The  mechanism 
consists  of  a  loading  button  glued  onto  a  rubber  diaphragm 
which  is  mounted  on  top  of  a  small  gas  chamber.  The 
mechanism  is  strategically  located  under  the  slider-spring 
assembly  such  that  when  the  small  gas  chamber  is 
pressurized,  the  diaphragm  expands  causing  the  loading 
button  to  protrude  and  to  push  against  the  slider  at  the 
pivot  point.  The  force  applied  to  the  slider  is  determine 
by  measuring  the  strain  on  the  cantilever  beam  to  which 
the  slider  assembly  is  attached.  The  beam  or  arm  is  made 
of  Delrin  material1  and  is  instrumented  with  four  strain 
gages.  Figure  11-10  shows  the  entire  arm  assembly.  The 
arm  is  calibrated  with  dead  weights  and  the  strain  rate  of 
change  per  gram  was  found  to  be  10.5  pm/m/gram  for  loads 
between  10  and  30  grams.  Figure  11-11  shows  the  typical 
calibration  curve  for  the  arm  at  20°C.  It  has  been  found 
that  the  calibration  curve  is  sensitive  to  ambient 
temperature.  Since  the  experiments  are  performed  in  an 
air-conditioned  laboratory,  close  control  of  the  ambient 
temperature  was  maintained. 

A  constant  temperature  hot-wire  anemometer  is  used  as 
a  gas  analyzer  to  measure  the  helium  gas  concentration 
inside  the  test  chamber.  Helium  was  selected  for  several 
reasons.  One  is  that  the  mean  free  path  of  helium  is  much 
larger  than  air  (0.188  inn  vs  0.069  wm) .  Another  reason  is 
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that  the  thermal  conductivity  of  helium  is  significantly 
different  (about  eight  times  higher)  than  air.  Thus  when 
helium  is  introduced  into  the  test  chamber ,  the  hot-wire 
anemometer  can  easily  detect  its  presence.  The  hot-wire 
probe  is  calibrated  inside  a  vacuum  bell  jar.  The  voltage 
response  of  the  probe  versus  the  helium  concentration  is 
plotted  in  Figure  11-12.  The  mean  free  path  curve  on  the 
figure  is  obtained  using  the  theoretical  formula  for 
mixtures  given  in  Reference  [293.  Therefore,  once  the 
probe  voltage  response  is  measured,  the  helium 
concentration  and  the  mean  free  path  of  the  helium/air 
mixture  are  easily  determined  using  the  calibration  curve. 
Viscosity  of  the  gas  mixture  can  also  be  calculated  using 
the  theoretical  mixture  formula  given  in  Reference  [293. 
It  should  be  noted  that  the  highest  purity  research  grade 
helium  is  used  in  the  experiments.  The  helium  is  also 
filtered  before  it  is  force  fed  into  the  test  chamber. 
The  hot  wire  probe  is  mounted  in  the  upper  chamber  and  is 
shielded  from  the  gas  motion  inside  the  chamber  caused  by 
the  rotating  disk.  Probe  responses  obtained  when  the  disk 
is  still  and  when  it  is  rotating  showed  no  difference. 
For  all  the  experiments,  the  test  chamber  was  filled  with 
lOOt  helium. 

The  entire  test  chamber  is  located  in  front  of  the 
laminar  flow  bench  in  order  to  prevent  particulate 
contamination  from  the  room  when  the  unit  is  opened.  In 
addition,  before  any  experiments  are  performed,  the  upper 
chamber  is  purged  with  clean  air  which  has  been  passed 
through  a  "clean  air"  filter. 


5.2  Optical  Bench 

The  optical  bench  used  to  measure  the  bearing 
clearances  is  separated  from  the  test  chamber.  The 
optical  bench  consists  of  two  monochromatic  light  sources 
of  different  wavelength,  a  beam  splitter,  3  mirrors,  beam 
expander,  focusing  lens,  and  an  imaging  system.  A 
schematic  of  the  interferometric  set-up  is  shown  in  Figure 
11-8.  Two  CW  lasers  are  used  to  generate  the  needed 
interference  fringe  patterns.  The  lasers  include  a 
Helium-Neon  laser  with  a  wavelength  of  632.8  nm  and  a 
Helium-Cadmium  laser  with  a  wavelength  of  441.6  nm.  Since 
the  final  approach  path  taken  by  the  two  laser  beams  must 
be  the  same,  and  since  the  lasers  cannot  be  placed 
collinear  with  each  other,  the  beams  are  made  collinear 
via  a  mirror/beam  splitter  set-up  as  shown  in  Figure  11-8. 
Once  the  beams  are  collinear,  they  are  expanded  to  about 
1.25  cm  in  diameter  via  a  beam  expander.  Due  to  the  path 
length  difference  between  the  beam  reflected  from  slider 
surface  and  that  reflected  from  the  bottom  surface  of  the 
glass  disk,  distinct  optical  interference  pattern  is 
observed  for  each  of  the  two  lasers.  These  fringe 
patterns  are  then  enlarged  ten-times  and  focused  onto  a 
Polaroid  film  back  where  black  and  white  photographs  are 
taken  for  each  laser.  Typical  photographs  of  the  fringes 
are  shown  in  Figure  11-13*  With  only  minor  adjustments, 
the  present  optical  arrangement  can  also  be  used  with 
white  light  (similar  to  Lin's  work  [22]). 


6.  EXPERIMENTAL  MEASUREMENTS 
6.1  Slider  Dimensions 

The  slider  bearings  used  in  this  experimental  study 
are  the  IBM  "Winchester"  type  flying  heads  found  in 
magnetic  disk  memory  packs.  Figure  1 1 —3  shows  a  sketch  of 
the  slider.  The  slider  is  made  of  sintered  ferrite  with 
its  two  skates  lapped  to  a  very  smooth  surface  finish. 

The  physical  dimensions  of  the  slider  are  obtained  by 
two  measuring  techniques.  An  optical  (machinist) 
microscope  with  10X  magnification  is  used  to  measure  the 
overall  dimensions  of  the  slider  head  and  its  skates. 
These  dimensions  include  the  length  and  width  of  the 
slider  head,  the  length  and  width  of  each  rail,  and  the 
projected  ramp  length  of  each  rail.  In  addition,  the 
microscope  is  also  used  to  measure  the  location  of  the 
load  point  (pivot  point).  Table  11-1  shows  the  dimensions 
for  all  the  various  sliders  used. 

Interferometric  techniques  are  used  to  measure  the 
taper  (ramp)  height  and  surface  crown  heights  because  they 
are  too  small  to  be  measured  by  the  microscope.  Both 
quantities  are  obtained  after  examining  the  ..interference 
fringe  photographs  taken  during  actual  experiments.  The 
ramp  height  is  determined  by  counting  the  number  of  dark 
fringes  appearing  across  the  ramp.  Like  all  the  other 
dimensions,  the  ramp  height  varies  from  slider  to  slider 
as  shown  in  Table  11-1.  Using  interferometry,  the  ramp 
height  can  be  determined  to  within  one-half  the  laser 
wavelength;  that  is,  for  He-Cd  laser,  the  error  is  only 
about  0.22  pm  out  of  9.0  pm. 


On  the  other  hand,  the  crown  height  on  the  slider 
skate  is  slightly  more  difficult  to  determine  from  the 
Interference  fringe  photographs.  Several  assumption  are 
made  before  the  crown  can  be  determined.  One  is  that  the 
crown,  when  present,  is  parabolic  in  shape  with  its  vertex 
located  at  the  center  of  the  "land"  section  of  the  skate. 
The  crown  may  be  either  concave  or  convex  relative  to  the 
flat  slider  and  is  easily  determined  by  examining  the 
manner  in  which  the  fringes  are  distributed  along  the 
skate.  Figure  11-15  illustrates  the  differences  among  the 
three  surface  contours  -  concave,  flat,  and  convex.  As 
illustrated,  for  a  concave  crown,  the  separation  between 
consecutive  fringes  becomes  larger  as  the  distance  from 
the  .trailing  edge  increases;  for  a  "flat"  surface,  the 
separation  remains  constant;  and  for  a  convex  crown,  the 
separation  becomes  smaller.  Another  assumption  made  is 
that  the  crowns  on  the  two  skates  do  not  necessarily  have 
to  be  the  same  in  magnitude  or  concavity.  The  last  column 
in  Table  11-1  shows  how  the  two  crowns  can  vary  on  each 
slider  head.  Furthermore,  it  is  also  assumed  that  the 
crown  does  not  change  in  shape  or  in  magnitude  within  the 
load  range  tested  (from  10  grams  to  30  grams).  Finally, 
it  is  assumed  that  there  is  no  crown  across  the  width  of 
the  skate. 

In  order  to  reduce  the  error  in  determining  the  crown 
height,  the  procedure  is  computerized.  In  a  previous 
study  [33],  the  crown  height  was  determined  manually  by 
fitting  an  assumed  crown  profile  through  data  points  from 
only  one  experiment.  Crown  shape  and  height  were  also 
assumed  to  be  the  same  on  both  skates.  However,  in  the 


present  study,  the  crown  concavity  and  height  are 
determined  based  on  all  the  experimental  data  for  the 
particular  slider  (usually  from  over  20  experiments). 
Since  the  crown  determination  is  directly  related  to  the 
determination  of  the  flying  profile,  the  exact  crowns  are 
determined  only  when  the  flying  profile  is  determined  for 
each  experiment.  The  exact  procedure  is  discussed  in  the 
following  section. 

6.2  Flying  Profile 

With  the  present  experimental  measuring  technique, 
the  actual  flying  profile  of  the  slider  can  be  determine 
with  great  accuracy.  The  procedure  used  to  determine  the 
profile  is  computerized  in  order  to  minimize  the  error  in 
data  reduction.  The  Method  of  Least  Squares  is  used  to 
determine  the  trailing  edge  clearance,  the  pitch  angle, 
and  the  roll  angle  for  a  given  slider  and  crown  profile. 
The  data  reduction  procedure  is  outlined  in  Appendix  C. 

Different  combination  of  the  crown  profiles  for  the 
two  skates  are  used  as  input  for  the  computer  fitting 
routine  and  the  combination  that  best  fits  the 
experimental  data  for  all  the  runs  (that  is,  gives  the 
least  (squares)  error  in  fitting)  determines  the  crown 
profiles  on  the  slider.  Moreover,  once  the  best  crown 
profile  combination  is  established,  the  "experimental" 
flying  profile  is  also  determined  for  each  experiment. 
This  new  data  reduction  procedure  is  more  accurate  than 
the  previous  method  used  [333  because  the  procedure  uses 
all  available  experimental  data.  Furthermore,  the  flying 
profile  determined  as  well  as  the  crown  are  mathematically 
the  best  for  the  given  data. 
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It.  has  been  determined  that  an  inaccurate  crown 
determination  will  lead  to  an  inaccurate  flying  profile 
determination.  With  this  new  procedure,  a  ±  0.025  ym 
accuracy  or  better  can  be  achieved  in  determining  the 
trailing  edge  clearance.  With  the  slider  bearing  expected 
to  operate  at  less  than  0.125  ym  minimum  film  thickness, 
an  error  of  more  than  0.025  ym  can  lead  to  significant 
errors  in  the  numerical  solution.  Previously,  the  best 

e 

accuracy  attained  was  no  better  than  ±  0.050  ym. 


7.  EXPERIMENTAL  RESULTS  AND  COMPARISON  TO  THEORY 

In  this  investigation  the  performance  of  Winchester 
type  heads  of  different  rail  width  -  0.38  mm,  0.51  mm,  and 
0.89  mm  -  were  studied.  The  physical  dimensions  of  the 
these  sliders  are  summarized  in  Table  11-1.  The  three 
types  of  slider  were  tested  in  both  normal  ambient 
conditions  (i.e.  air  at  room  conditions)  and  in  a  pure 
helium  environment.  Whenever  possible,  the  same  slider  is 
used  in  both  gas  media.  However,  this  was  not  possible  in 
most  cases  because  the  slider  has  the  tendency  to  "crash" 
at  high  loads  and  low  clearances  due  to  dust 
contamination.  Only  the  wide  head  (0.89  mm  wide)  did  not 
crash.  Crash  is  defined  as  the  slider  surface  being 
scratched  by  dust  trapped  in  the  gas  film.  The  slider 
will  still  "fly"  but  no  data  can  be  taken  because  the 

scratched  slider  surface  makes  the  interference  fringes 
impossible  to  identify.  Figure  11-16  shows  what  happens 
to  the  fringes  when  the  slider  "crashes".  (Note  the 

streaking  along  the  length  of  shoe  B. ) 

Each  type  of  slider  was  tested  at  various  velocities 
between  17  to  52  m/sec  with  loads  between  8  to  30  grams. 
The  clearance  profile  for  each  operating  condition  is 
determined  using  the  method  of  least  squares  for  fitting 
the  experimental  data  as  describe  in  Appendix  C.  The 
experimental  results  are  then  compared  with  numerical 

solutions  of  the  modified  Reynolds  equation. 

The  theoretical  prediction  are  based  on  actual 
physical  dimensions  of  the  slider  and  not  on  the  nominal 
design  values.  Gas  properties  of  air  and  helium  are 

assumed  to  be  constant  throughout  the  experiments;  that  is 
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air  and  helium  viscosities  are  taken  to  be  18.27  Pa-s  and 
19.50  Pa-s,  respectively,  while  the  air  and  helium  mean 
free  paths  are  taken  to  be  0.069  pm  and  0.188  pm, 
respectively.  The  ambient  pressure  needed  for  calculating 
the  theoretical  load  is  taken  to  be  101.4  kPa. 


There  are  two  methods  by  which  the  results  can  be 
compared.  The  first  method  is  to  compare  the  accuracy  of 
the  predicted  trailing  edge  clearance  and  experimentally 
measured  clearance  for  a  given  load.  With  this  approach, 
the  pivot  position  of  the  slider  is  assumed  to  be  known. 
Armed  with  this  assumption  and  a  starting  reference 
clearance  (usually  it  is  the  trailing  edge  clearance),  the 
computer  program  proceeds  to  solve  for  the  flying  profile 
whose  center  of  pressure  coincides  with  the  pivot  point. 
The  program  also  calculates  the  load  associated  with  the 
profile.  In  theory,  if  the  load  does  not  correspond  to 
the  given  load,  the  program  will  change  the  reference 
clearance  and  repeat  the  calculations  until  both  the  load 
and  pivot  point  requirements  are  satisfied.  In  practice, 
this  type  of  iterative  scheme  is  not  used  because  enormous 
an  amount  of  computational  time  is  needed  to  solve  just 
one  case.  Instead,  the  program  will  generate  a  load 
versus  trailing  edge  clearance  curve  for  fixed  pivot 
point.  Thus,  to  determine  the  theoretical  clearance  for  a 
given  load,  one  only  has  to  go  to  the  curve.  Associated 
with  the  load/spacing  curve  is  a  pitch/spacing  curve.  In 
this  investigation,  the  measured  and  the  predicted  pitch 
angles  is  also  studied. 

The  second  method  is  to  compare  the  predicted  load 
and  the  predicted  center  of  pressure  for  given  flying 
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profile  (i.e.  specify  the  trailing  edge  clearance  and  the 
pitch  angle) . 

The  experimental  curves  for  the  load/spacing  and 
pitch  angle/spacing  characteristics  of  the  three  different 
width  sliders  in  air  and  helium  are  presented  in  Figures 
11-17  to  11-32,  Figures  11-33  to  11-50,  and  Figures  11-51 
to  11-62.  Also  plotted  on  these  figures  are  the  numerical 
solutions  of  the  modified  Reynolds  equation.  All  the 
theoretical  curves  are  for  a  surface  accommodation 
coefficient  of  0.89,  unless  otherwise  indicated.  No  roll 
angle/spacing  characteristics  were  plotted  because  very 
small  roll  angles  were  encountered  even  when  the  two  shoes 
on  the  sliders  have  different  surface  curvatures. 
Furthermore,  even  if  there  was  a  small  roll  angle,  the 
predicted  load  and  pitch  did  not  differ  by  more  than  5% 
when  the  roll  angle  was  taken  into  account.  Since  three 
different  width  sliders  were  tested,  in  order  to  avoid 
confusion,  the  sliders  are  classified  into  three 
categories  according  to  their  widths  -  narrow,  standard, 
and  wide. 

7.1  Narrow  Sliders 

Narrow  sliders  were  tested  in  air  and  helium  at 
various  disk  velocities.  The  slider  used  in  the  air 
experiments  was  not  used  in  the  helium  experiments  because 
it  had  crashed  at  the  end  of  the  air  runs.  Another  slider 
of  the  same  nominal  dimensions  was  used  in  helium.  Figure 
11-17  and  11-18  show  the  typical  experimental  load/spacing 
relationship  for  the  narrow  sliders  in  air  and  helium, 
respectively,  at  the  same  disk  velocity.  For  each 
load/spacing  plot,  there  is  a  corresponding  pitch 


angle/spacing  plot.  Figures  11-19  and  11-20  show  the 
pitch/ spacing  characteristics  for  the  two  cases  presented 
in  Figures  11-17  and  11-18,  respectively.  All  theoretical 
curves  are  based  on  surface  accommodation  coefficient  of 
0.89  for  both  air  and  helium  cases,  unless  otherwise 
indicated . 

Several  important  observations  can  be  made  from 
examining  Figures  11-17  and  11-18.  To  illustrate  the 
effects  of  molecular  slip,  theoretical  curves  based  on 
continuum  theory  and  slip  flow  theory  are  plotted  on  the 
two  figures.  As  the  figures  clearly  indicate,  the  main 
effect  of  slip  is  to  decrease  the  load  carrying  capacity 
of  the  slider  bearing  since  the  slip  theory  curve  always 
lies  below  the  continuum  theory  curves.  Furthermore,  it 
is  observed  that  the  effect  of  slip  is  greater  in  helium 
than  air  demonstrated  by  the  fact  that  the  separation 
between  slip  theory  curve  and  continuum  theory  curve  is 
greater  in  helium  than  in  air.  This  is  attributed  to  the 
fact  that  the  mean  free  path  of  helium  is  much  larger  than 
that  of  air.  In  addition,  the  two  figures  also  show  that 
even  though  the  conventional  bearing  number  (indicated  by 
the  upper  scale  on  the  figure)  is  of  the  order  10^  to  101*, 
the  effects  of  slip  and  side  leakage  cannot  be  neglected. 
If  the  high  bearing  number  effect  had  dominated,  the 
slip-flow  theory  curves  and  the  continuum  theory  curves 
would  have  coincided,  as  was  the  case  reported  by  Tseng 
[35].  At  first  this  may  seem  to  be  a  surprising  result. 
However,  if  the  modified  bearing  number  as  defined  by 
equation  (11.24)  is  used,  it  is  obvious  that  the  bearing 
is  far  from  approaching  the  high  "bearing"  number  range. 
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In  fact,  for  the  narrow  slider,  the  modified  bearing 
number  is  only  on  the  order  of  101  to  10^. 

Figures  11-17  and  11-18  also  show  the  effect  of 
different  values  of  surface  accommodation  coefficient,  a  , 
on  the  theoretical  predictions.  Two  values  of  a  were 

used  -  1.0  and  0.89.  As  the  figures  show,  slip  theory 

with  a  =  0.89  gives  a  slightly  better  agreement  with 
experiment  for  both  air  and  helium.  The  surface 
accommodation  coefficient  for  air/glass  interface  was 
measured  to  be  0.89  [30],  and  Tseng  [35]  has  shown  that 
when  a  is  taken  as  0.89,  a  better  agreement  between 
theory  and  experiment  can  be  achieved  for  air.  However, 

there  is  no  published  data  on  o  for  helium/glass 

interface.  In  this  investigation,  different  values  of 
a  were  tried  and  it  was  found  that  when  o  is  equal  to 
0.89  a  good  agreement  between  experiment  and  theory  can  be 
achieved . 

Further  examination  of  Figures  11-17  and  11-18  shows 
that  as  the  clearance  is  decreased  by  increasing  the 
external  load,  the  agreement  between  experiment  and  theory 
does  not  deteriorate  even  though  the  ambient  Knudsen 
number  is  increased  from  0.1  to  0.38  for  air  and  from  0.6 
to  1.1-4  for  helium.  But,  it  should  be  noted  that  much 
closer  agreement  is  found  for  helium  experiments  than  for 
air  runs.  Moreover,  for  air,  Figure  11-17  also  shows  some 
departure  between  experiment  and  slip-flow  theory  at  lower 
spacings.  It  is  believed  that  helium  results  are  better 
than  air  because  helium  is  monatomic  while  air  is  a 
mixture  polyatomic  gases.  Thus  the  ideal  gas  assumption 
and  the  slip  flow  approximation  used  to  derive  the 


modified  Reynolds  equation  may  be  more  applicable  for 
helium  than  for  air,  especially  up  in  the  higher  Knudsen 
number  ranges. 

Examination  of  Figures  11*19  and  11-20  shows  that  the 
slip  flow  theory  with  o  =  0.89  is  also  able  to  predict 

with  accuracy  the  pitch  angle  of  the  slider  bearing 

operating  in  air  and  helium. 

Results  for  the  narrow  slider  at  other  operating 

conditions  are  presented  in  Figures  11-21  to  11-26  and 
Figures  11-29  to  11-32  for  air  and  helium,  respectively. 
It  is  clear  from  examining  these  figures  that  slip  flow 
theory  with  o  =  0.89  is  able  to  predict  with  accuracy  the 

bearing  performance  over  a  wide  range  of  bearing  numbers 
and  for  high  ambient  Knudsen  numbers. 

7.2  Standard  Sliders 

Experimental  results  for  the  standard  slider 
operating  at  various  velocities  in  air  and  helium  are 

presented  in  Figures  11-33  to  11-38  and  in  Figures  11-39 
to  11-50,  respectively.  Once  again,  slip  flow  theory  with 
o  s  0.89  seems  to  give  the  best  predictions  for  both  air 
and  helium  cases.  Furthermore,  it  is  found  that  slip  flow 
approximation  does  not  fail  even  when  the  ambient  Knudsen 
number  is  much  larger  than  1.0.  In  particular,  Figure 
11-39  shows  that  the  slip  flow  theory  is  able  to 
accurately  predict  the  load/spacing  characteristics  of  the 
slider  bearing  even  down  to  a  trailing  edge  clearance  of 
0.075  pm  with  the  ambient  Knudsen  number  of  2.5.  It 
should  be  noted  that  this  agreement  is  not  due  to  the  high 
bearing  number  effect  because,  once  again,  the  continuum 
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theory  curve  and  the  slip-flow  theory  curve  do  not 
coincide  as  shown  by  Figure  11-39. 

It  should  be  noted  that,  similar  to  the  narrow 
slider,  the  agreement  between  experiment  and  theory  is 
better  for  helium  cases  than  for  air.  Moreover,  the 
departure  between  experiment  and  theory  for  clearances 
below  0.25  ym  -  which  was  found  in  the  narrow  slider 
results  for  air  -  is  also  present  in  the  standard  slider 
air  results. 

Once  again,  Figures  11-36  to  11-38  and  Figures  11-45 
to  11-50  show  that  the  slip  flow  theory  with  as  0.89  can 
also  predict  with  accuracy  the  pitch  angle  of  the  standard 
slider  in  air  and  in  helium,  respectively. 

7.3  Wide  Sliders 

Experimental  results  for  the  wide  slider  operating  at 
various  velocities  in  air  and  in  helium  are  presented  in 
Figures  11-51  to  11-53  and  Figures  11-57  to  11-59, 

respectively.  The  same  slider  head  was  used  in  air  and  in 

helium  experiments.  Once  again,  the  slip  flow  theory  with 
o  =  0 . 89  seems  to  give  a  better  agreement  between 
experiment  and  theory  for  both  air  and  helium.  However, 
the  experimental  results  for  the  wide  head  are  more 

scattered  than  those  for  either  the  narrow  or  the  standard 
heads.  This  scatter  can  be  explained  by  the  fact  that  the 
surface  on  the  wide  head  is  very  irregular  and  the  surface 
profile  assumed  for  data  fitting  was  only  a  rough 
approximation  of  the  real  contour.  This  deficiency  is 
more  evident  when  comparing  the  measured  pitch  with 

predicted  pitch.  Figures  11-54  to  11-56  and  Figures  11-60 


to  11-62  show  the  inability  of  the  theory  to  predict  the 
slider  pitch  angle  for  both  air  and  helium.  It  must  be 
noted  that  this  failure  is  due  to  the  inability  to  account 
for  the  surface  irregularities  and  not  due  to  failure  of 
the  slip-flow  model  because  the  ambient  Knudsen  numbers 
are  all  below  0. 1 . 

Comparing  the  load/spacing  results  from  the  three 
different  width  sliders,  one  sees  that  the  wide  slider 
operates  at  much  higher  clearances  than  the  narrow  slider 
even  though  the  width  of  the  wide  bearing  is  only  twice 
that  of  the  narrow  bearing. 


no 


8.  DISCUSSION 

Results  presented  in  the  previous  chapter  clearly 
indicate  that  the  modified  Reynolds  equation  with  slip 
flow  approximation  can  predict  correctly  the  reduction  in 
load  capacity  of  finite  width  slider  bearings.  In 
addition,  the  slip-flow  theory  can  also  predict  accurately 
the  flying  pitch  angle.  The  agreement  between  experiment 
and  slip  flow  theory  is  excellent  for  low  ambient  Knudsen 
numbers  (Kg,  <  0.1),  moderate  Knudsen  numbers 
(0.1  <  K«  <  0.5),  and  high  Knudsen  numbers 
(0.5  <  Koo  <  2.5).  It  is  somewhat  misleading  to  use  the 
ambient  Knudsen  number  as  the  measure  for  molecular  slip 
effect  because  the  mean  free  path  of  the  gas  inside  the 
gas  film  is  smaller  than  the  ambient  value  due  to  the  fact 
that  the  local  pressure  is  higher  than  the  ambient 
pressure.  Thus,  the  local  Knudsen  number,  Kj ,  based  on 
local  values  of  mean  free  path  and  clearance  would  give 
more  insight  into  the  molecular  slip  effect.  Figure  11-63 
and  Figure  11-6M  are  typical  contour  plots  of  constant 
local  Knudsen  numbers  for  the  narrow  slider  in  air  and 
standard  slider  in  helium,  respectively.  As  the  figures 
clearly  show,  in  terms  of  local  Knudsen  numbers,  most  of 
the  bearing  is  still  dominated  by  relatively  high  local 
Knudsen  numbers.  Based  on  these  typical  plots  of  local 
Knudsen  numbers,  it  is  clear  that  the  effect  of  molecular 
slip  is  still  substantial  even  with  bearing  numbers  of 
6090  and  19308  for  the  two  cases  shown,  respectively. 

In  this  investigation,  the  effect  of  slip  has  been 
successfully  studied  without  becoming  concerned  with  high 
bearing  number  effect.  Had  the  bearing  number  been  large, 
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the  effect  of  slip  would  be  diminished.  Moreover,  there 
would  be  a  small  boundary  layer  at  the  side  edges  of  the 
slider.  To  show  that  high  bearing  number  limit  is  not 
reached  even  though  the  conventional  bearing  number  is 
large  for  most  of  the  experimental  runs,  pressure  profiles 
across  the  slider  bearing  are  plotted  for  different  points 
along  the  length  of  the  bearing.  Figure  11-65  and  11-66 
are  typical  plots  of  the  pressure  profiles  across  the 
bearing  for  the  narrow  slider  in  air  and  for  the  standard 
slider  in  helium,  respectively.  The  conventional  bearing 
number  for  Figure  11-65  is  6091  and  for  Figure  11-66  is 

19308.  It  is  clear  that  neither  of  the  pressure  profiles 

exhibit  any  sign  of  a  boundary  layer  at  the  side  edges 
because  there  are  no  sharp  pressure  gradients  near  the 
edges . 

In  this  investigation,  the  agreement  between 

experiment  and  slip-flow  theory  is  consistently  better  for 
helium  results  than  air  results.  It  is  believed  that  this 
phenomenon  is  caused  by  the  fact  that  helium  is  a 

monatomic  gas  while  air  is  a  mixture  of  polyatomic  gases. 
Thus  the  slip-flow  theory  derived  on  the  assumption  of 
ideal  gas  may  be  more  applicable  for  helium  than  for  air 
especially  in  the  higher  Knudsen  number  range,  i.e.  in 
transition  flow  regime.  It  was  noted  by  Bird  [1]  that  if 
a  gas  flow  is  in  local  thermodynamic  equilibrium,  the 
question  of  whether  it  is  a  simple  gas  or  gas  mixture  is 
of  no  consequence.  Bird  also  asserted  that  in  transition 
flow,  a  gas  mixture  may  be  initially  in  uniform 
composition,  but  species  separation  may  occur  due  to 
thermal  or  pressure  diffusion.  Whether  this  separation  is 


present  in  the  present  study  is  not  known.  A  sore 
comprehensive  analysis  of  this  phenomenon  is  needed  in 
order  to  fully  understand  why  helium  results  are  better 
than  air  results. 

From  examining  the  load/spacing  characteristic 
curves,  it  appears  that  slip  flow  theory  can  be  made  to 
agree  better  with  the  experimental  data  by  choosing  the 
surface  accommodation  coefficient,  a  ,  to  be  0.89  for  both 
air/glass  and  helium/glass  interfaces.  Since  the  surface 
accommodation  coefficient  is  directly  related  to  the 
surface  finish  of  the  solid  surface,  one  could,  perhaps, 
interpret  this  coefficient  as  a  crude  measure  of  surface 
roughness  with  amplitude  of  less  than  one  mean  free  path 
of  the  gas.  For  diffuse  surface  like  a  rough  surface,  the 
molecules  hitting  the  surface  are  fully  accommodated  and 
o  s  1.0.  On  the  other  hand,  if  the  surface  is  not  totally 
diffuse  like  a  shinny  surface,  the  molecules  hitting  the 
surface  are  reflected  specularly  and  diffusely  and 
a  <  1.0.  Therefore,  if  the  predictions  are  plotted  for 
two  values  of  o  ,  the  curve  for  lower  value  of  o  would 
lie  below  the  higher  value  of  o  as  shown  in  Figure  11-17. 
Thus  in  terms  of  "roughness",  the  smoother-finish  (shinny) 
surface  corresponding  to  the  lower  value  of  o  would  carry 
less  load  than  the  rough  surface  corresponding  to  the 
higher  value  of  o  .  This  result  is  similar  to  that 
obtained  through  formal  treatment  of  surface  roughness. 
It  should  be  noted  that  this  is  only  an  interpretation  of 
surface  accommodation  coefficient  in  terms  of  surface 
finish . 

The  agreement  observed  between  experiment  and  theory 
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also  suggests  that  the  numerical  scheme  used  to  solve  the 
Reynolds  equation  is  both  reliable  and  accurate.  In  this 
investigation,  several  different  numerical  algorithms  were 
used  to  solve  the  modified  Reynolds  equation;  for  example, 
integral  discretization  method  with  variable  grid  and 
exponential  approximation  with  constant  and  variable 
grids.  The  variable  grid  used  in  the  numerical  algorithms 
is  dependent  on  the  operating  conditions  and  is  describe 
in  Appendix  E.  All  the  algorithms  used  the  Column  Method 
of  solution  for  solving  the  system  of  equations  resulting 
from  the  finite  difference  approximation.  The  purpose  of 
this  exercise  was  to  determine  the  most  efficient 
algorithm,  in  terms  of  solution  accuracy  and  computational 
time,  for  solving  the  modified  Reynolds  equation.  Table 
11-67  shows  the  comparison  of  the  three  algorithms  used. 
It  is  quite  clear  from  Table  11-67  that  the  slowest  but 
very  accurate  algorithm  is  the  integral  discretization 
method  with  variable  grid  while  the  fastest  but  not  very 
accurate  algorithm  is  the  exponential  approximation  with 
constant  grid.  But  the  exponential  approximation  with 
variable  grid  is  both  fast  and  accurate.  It  is  60  percent 
faster  in  computational  time  and  is  within  2%  of  the 
predicted  results  (load  and  pitch  angle)  when  compared 
with  the  integral  discretization  method  with  variable 
grid.  It  should  be  noted  that  variable  grid  is  used  in 
the  exponential  approximation  method  to  improve  solution 
accuracy  and  not  for  numerical  stability. 

Numerical  algorithms  using  the  Alternating 
Directional  Implicit  (ADI)  scheme  to  solve  the  system  of 
equations  resulting  from  two  discretization  methods  - 


integral  discretization  and  exponential  approximation  - 

were  also  examined  because  in  general  algorithms  using  ADI 

technique  is  usually  faster  than  other  schemes  using 

matrix  inversions  like  the  Column  Method.  To  add 

stability  to  the  ADI  scheme,  the  time  dependent  term  in 

the  Reynolds  equation,  MUilL  ,  is  retained  and  the  problem 

3T 

is  solved  as  a  transient  one.  The  Peaceman-Rachford  ADI 
method  [28]  was  used  and  is  outlined  the  Appendix  D.  In 
this  study,  Successive  Line  Relaxation  method  with 
Alternating  Direction  was  also  used  and  is  also  discussed 
in  Appendix  D. 

For  algorithms  using  integral  discretization 
technique,  it  was  found  that  the  ADI  solution  scheme  was 
much  slower  than  the  Column  Method  because  the  size  of  the 
time  step  was  limited  by  the  smallest  grid  spacing  used. 
It  was  found  that  for  each  problem  there  was  an  optimal 
size  for  the  time  step  with  which  a  stable  steady  state 
solution  can  be  obtained  with  the  minimum  amount  of 
iterations.  Therefore,  if  the  time  step  was  too  large, 
the  algorithm  would  become  unstable  and  no  solution  could 
be  obtained;  on  the  other  hand,  if  the  time  step  was  too 
small,  the  algorithm  would  take  forever  to  reach  the 
steady  state  solution.  However,  for  the  test  case 
studied,  even  with  the  optimal  time  step,  the  ADI  scheme 
was  still  slower  than  the  Column  Method. 

It  has  also  been  found  that  algorithm  using  the 
exponential  approximation  technique  with  ADI  was  extremely 
unstable  for  both  constant  and  variable  grids.  The 
algorithm  was  found  to  be  stable  only  when  the  time  step 
was  at  least  two  orders  of  magnitude  smaller  than  that 


needed  for  the  algorithm  using  integral  discretization 
with  ADI.  Once  again,  there  seemed  to  be  an  optimal  time 
step.  However,  if  Alternating  Directional  Successive  Line 
Relaxation  method  was  used  with  the  exponential 
approximation  technique,  the  algorithm  was  found  to  be 
unconditionally  stable  for  any  size  time  step  with  either 
constant  or  variable  grid  spacing.  However,  even  for  an 
infinite  time  step,  the  algorithm  was  found  to  be  still 
several  seconds  slower  than  same  algorithm  with  Column 
Method . 
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9.  GENERAL  REMARKS  AND  RECOMMENDATIONS 

The  validity  of  the  modified  Reynolds  equation  with 
slip  flow  approximation  has  been  verified  over  the  entire 
range  of  Knudsen  numbers  between  0  and  1.0. 

Nevertheless,  the  author  believes  that  additional 
research  is  need  in  the  following  areas: 


*  Determine  the  reasons  why  helium  results  are 
better  than  air 


*  Develop  faster  numerical  algorithms  for 
solving  the  modified  Reynolds  equation 


*  Determine  the  extent  of  surface  roughness 
effects  in  narrow  bearing  systems 


10.  NOMENCLATURE 

h  Film  thickness 

ht  Trailing  edge  film  thickness 

hm£n  Minimum  film  thickness 

hr  Ramp  height 

hc  Crown  height  of  slider  shoe 

H  Dimensionless  film  thickness  (  =  h/ h,,, ^ n  ) 

Ht  Dimensionless  trailing  edge  film  thickness 

(  =  bt/hmin  ) 


K 

K 

K 

1 

1 

1 


1 


r 

piv 


P 

Pa 

P 

Q 


w 

U 


X 

y 

t 

x 


Y 


Knudsen  number  (  s  X /h  ) 

Ambient  Knudsen  number  (  =  x  a/ in  ) 

Local  knudsen  number  (  =  X  j/h  ) 

Length  of  slider 
Ramp  length  (projected) 

Pivot  point  location  (measured  from  the 
trailing  edge) 

Pressure 

Ambient  pressure 

Dimensionless  pressure  (  s  p/pa  ) 
Linearized  variable  (  s  P^H^  ) 

Width  of  slider  shoes 
Linear  velocity  of  disk 

Components  of  velocity  in  x,y  directions, 
respectively 

Spatial  coordinate  in  length  direction 
Spatial  coordinate  in  width  direction 
Time 

Dimensionless  spatial  coordinate  along  the 
length  (  =  x/1  ) 

Dimensionless  spatial  coordinate  along  the 
width  (  s  y/w  ) 


W 


V 


Load 

Viscosity  of  gas 
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r  !c 


4 

;  c 

•L 


Subscripts: 


a  ,  - 


ith  grid  spacing  in  X  direction 

jth  grid  spacing  in  Y  direction 

Wavelength  of  monochromatic  light  source 
or  molecular  mean  free  path  of  the  gas 

Bearing  number  based  on  trailing  edge 
clearance,  (  «  6uUl/paht2  ) 

Surface  accommodation  coefficient 

Numerical  constant  (usually  set  equal  to 


ambient  conditions 

x,y  index  of  grid  point,  respectively 
trailing  edge 
ramp  region 

local  conditions  or  values 
minimum  value 
pivot  point 
crown 

x,y  components  of  vector  quantity 


SLIDERS  TESTED 


(b)  END  VIEW 
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mole  fraction  of  air 


FIGURE  11-12:  HOT  WIRE  PROBE  RESPONSE  CURVE  FOR 
HELIUM/AIR  MIXTRUE 


MIXTURE  MEAN  FREE  PATH,  pm 


PA3  SLIDER  in  HELIUM  (SHOE  A) 


FIGURE  11-14:  DATA  PROCESSING  AND  CROWN  DETERMINATION 


PHENOMENON  -  PA 3  SLIDER  IN  AIR  -  U  -  36.07  m/s^c 


CLEARANC 


PITCH  ANGLE  VS  TRAILING  EDGE  CLEARANCE,  U  *  36.13  m/sec 


16618 


FIGURE  H -21:  LOAD  VS  TRAILING  EDGE  CLEARANCE ,  I  =  44.88  m/sec 


figure  lt-24:  LOAD  VS  TRAILING  EfGE  CLEAR ANC 


7661  4415  1962 


MGURE  11-29:  LOAD  VS  TRAILING  EDGE  CLEARANCE,  U  =  44.88  m/sec 


FIGURE  t 1 -30:  LOAD  VS  TRAILING  EDGE  CLEARANCE ,  U  =  51.90  m/sec 


PITCH  ANGLE  VS  TRAILING  EDGE  CLEARANCE,  U  =  44.88  m/sec 


VS  TRAILING  EDGE  CLEARANCE,  U  =  44.88  m/sec 


rEARANCH ,  U  =  51.90  tn/sec 


PITCH  ANGLE  VS  TRAILING  EDGE  CLEARANCE,  U  =  44.88  m/sec 


PITCH  ANCLE  VS  TRAILING  EDGE  CLEARANCE,  U  =  51.90  m/sec 


FIGURE  11-40:  load  vs  trailing  edge  CLEARANv.:-  u  =  2U.59  m/sec 


FIGURE  11-41:  LOAD  VS  TRAILING  EDGE  CLEARANCE,  U  =  30.53  m/sec 


TRAILIN'".  i:nCE  CLEARANCE,  ll  =  36.71  m/sec 


LOAD  VS  TRAILING  EDGE  CLEARANCE/  U  =  m/sec 


2058*  5146  2287  1286  821  572 


TRAILING  EDGE  CLEARANCE 


12170  30*2  1352 _ 761 _ 087 _ 338 


PITCH  ANGLE  VS  TRAILING  EDGE  CLEARANCE,  U  =  30.54  m/sec 


»632 _ 3658 _ 1626 _ 914 _ 585 _ 406 


PITCH  ANGLE  VS  TRAILING  EDGE  CLEARANCE,  U  *  36.71  m/sec 


7764  «<MM 


PITCH  ANGLE  VS  TRAILING  EDGE  CLEARANCE,  U  =  44.56  m/sec 


r  IGURE  11-50:  PITCH  ANGLE  VS  TRAILING  EDGE  CLEARANCE,  0  =  51.63  m/sec 


I.OAP  VS  TRAILING  EDGE  CLEARANCE,  U  =  30.13  m/soc 


r  IGURE  11-53:  LOAD  VS  TRAILIN'';  EDGL  CLKARANCL ,  U  =  r,1.90  m/sec 


259 _ 205 _ 166 _ 137 


PITCH  ANGLE  VS  TRAILING  EDGE  CLEARANCE,  U  =  44.88  m/sec 


PITCH  ANGLE  VS  TRAILING  EDGE  CLEARANCE 


FIGURE  11-57:  LOAD  VS  TRAILING  EDGE  CLEARANCE,  U  *  35.35  m/sec 


247  200  166  139 


PITCH  ANGLE  VS  TRAILING  EDGE  CLEARANCE,  U  *  50.81  m/sec 
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LEADING  EDGE 


0.8t 


K  »  0.10 


0.6L 


FIGURE  11-64:  CONTOUR  PLOT  OF  LOCAL  KNUDSEN  NUMBER 
UNDER  A  SLIDER  PAD  (  Y2  SLIDER  at 
ht  -  0.075  um,  a  -  19308,  •  2.47) 


PA3  SLIDER:  (SHOE  A) 
(NARROW) 

ht  *  0.1 905  urn 

6p  *  38.92  urad 
A  *  6091 
K  *  0.35 
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DERIVATION  OF  THE  MODIFIED  REYNOLDS  EQUATION 


The  modified  Reynolds  equation  as  given  by  equation 
C 3 - 3 >  is  derived  in  this  appendix.  The  equation  was  first 
developed  by  Burgdorfer  [53.  In  this  investigation ,  the 
equation  is  slightly  different  from  Burgdorfer's  to  give  a 
better  insight  into  the  effect  of  side  leakage.  The 
equation  below  are  written  with  reference  to  the  cartesian 
coordinate  system  defined  in  Figure  11-2. 

The  appropriate  Navier-Stokes  equations  and 
continuity  equation  for  the  gas  bearing  problem  are: 


Momentum : 


(ii.i) 

(11.2) 

H-o 

(11.3) 

Continuity: 

,£ud*]  +  •  0 

(11. A) 

The  boundary  conditions  are: 

u 

(s-0)  -  U  +  X 

a*  |z«o 

(11.5) 

o 
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( 


( 


C 


c 


c 


t 


u  (z«h)  -  -X 

a2lz«h 

(11.6) 

v  (z-0)  -  X 

I  z«=0 

(11.7) 

v  (z-h)  -  -X  fl|2.h 

(11.8) 

From  equation  (11.3),  it  follows  that  the  pressure  p 
is  independent  of  z.  Integrating  equation  (11.1)  twice 

with  respect  to  z  and  using  the  boundary  conditions  (11.5) 

and  (11.6)  gives: 

u '  hz  -  *2]  +  u[ 

,  z  +  X  1 

Similarly,  solving  equation  (11.2) 

(11.7)  and  (11.8)  gives: 

with  boundary  condition 

v  -  -  jp  |?[hX  +  hz  -  *2] 

(11.10) 

Substituting  equations  (11.9)  and  (11.10)  into  equation 
(11 .  M)  gives: 


(11.11) 

Because  the  viscosity  is  nearly  Independent  of 
pressure  in  the  range  considered  and  since  the  teaperature 
variations  across  the  gas  film  are  saall,  the  viscosity 
oan  be  assuaed  to  be  constant.  Moreover,  an  isotheraal 
relationship  can  be  assuaed,  l.e. 


(11.12) 


With  the  above  assumptions,  equation  (11.11)  becomes 


£[ph3H(l  +  TT>]  +  +  TT>]  -  <»4&<ph>] 

(11.13) 

Equations  (11.11)  and  (11.13)  are  modified  forms 
Reynolds  equation  for  the  slip  flow  region  (0  <  A /h  <  1 ) . 
The  boundary  condition  on  equation  (11.11)  or  (11.13)  is: 


P(x,y)  =  pa 

on  all 

boundaries  (11.14) 

To 

write  equation 

(11.13) 

in  dimensionless  form, 

define 

the  following: 

P  = 

P/Pa 

(11.15) 

X  = 

x/1 

(11.16) 

Y  = 

y/w 

(11. 17) 

H  = 

h/hmin 

(11.18) 

For 

constant  temperature, 

A 

Pa 

*a 

P 

(11.19) 

With  equations  (11.15)  to  (11.19),  equation  (11.13) 
becomes : 


(t)  Mph3|*(1  +  tf f}]  +  Mph3|?(1  +  i?w)]  "  a(t^(ph) 


(11.20) 


6jjU1 


where 


(11.22) 


In  equation  (11.20),  the  slider  width  to  length  ratio 
appears  as  the  result  of  the  non-dimensional ization 
process.  With  this  ratio  present,  the  effect  of  side 
leakage  is  more  evident.  As  the  equation  clearly  shows, 
the  side  flow  term  (i.e.  the  Y-derivative  term)  may  not  be 
neglected  if  the  width  to  length  ratio  is  small  since  the 
other  terms  in  the  equation  are  multiplied  by  the  factor 
(w/l)^.  Typically,  the  (w/1)  ratio  for  the  present  narrow 
bearing  is  about  0.1. 


Equation  (11.20)  can  be  written  in  the  following 


form: 


(?)2£H&'  ♦  £>]  ♦  *  &>]  -  ^ 


(11.23) 


where  the  Modified  Bearing  Number  is  define  as  - 


A  -  A 


(11.24) 


APPENDIX  B 


EXPONENTIAL  APPROXIMATION  OF  REYNOLDS  EQUATION 


In  this  appendix,  the  exponential  approximation  of 
the  Modified  Reynolds  equation  is  outlined.  This 
technique  was  originally  developed  by  H.G.  Elrod  [12]. 
The  author  is  very  grateful  to  Professor  Elrod  for  sharing 
this  algorithm. 

The  steady-state  modified  Reynolds  equation  can  often 
be  written  in  terms  of  the  mass  fluxes: 


w'"1*1  +  »Y(,V  ‘  ° 


(11.25) 


where 


-  PH3(1 


*  TTT)|?  +  AxHP 


(11.26) 


my  -  -  PH3(1  +  +  AyHP 


(11.27) 


with  the  variables  nondimensionalized  as  follows: 


X 


Now  consider  the  mass  flux  in  the  X  direction  as 
given  by  equation  (11.26),  which  is  is  a  non-linear  first 
order  differential  equation  in  P.  Assuming  this  mass  flux 
(mx)  is  constant  between  two  points,  say  and  Xi+1,  and 


also  assuming  the  coefficients  in  equation  (11.26)  are 
constant  with  some  average  value  between  the  two  points, 
Equation  (11.26)  can  be  rewritten  as: 


+  C 


2 


where 

C1  -  PH3  (1  +  ^0 


( 1 1 . 28 ) 

(11.29) 


C 


2 


(11.30) 


(The  bar  denotes  some  averaged  value  between  X^  and  Xi+1) 

The  boundary  condition  on  equation  (11.26)  is: 

P  :  Pt  at  X  =  XL  (11.3D 

The  solution  for  equation  (11.26)  with  the  above 
boundary  condition  (11. 31)  is: 


or  solving  for  mx ,  equation  (11.32)  becomes  - 


(11.32) 


+  c2p, 


(11.33) 


Therefore,  the  mass  flux  between  X^  and  X^+1  is  given  by 


where  Ci  and  C2  are  given  by  equations  (11.29)  and 
(11.30),  respectively. 

Similar  expression  can  be  obtained  for  the  mass  flux 
in  the  Y  direction,  i.e.  between  Yj  and  Yj+1: 


c 


o 


c 


(11.35) 


(11.36) 


C 


4 


(11.37) 


(The  bar  denotes  some  average  value  between  Yj  and  Yj+^) 

With  mx  and  <hy  known,  standard  finite  difference 
approximation  can  be  applied  to  equation  (11.25)  for  each 
grid  point  and  the  problem  can  thus  be  solved  with  either 
Column  Method  or  Alternating  Directional  Implicit  (ADI) 
scheme  or  Alternating  Directional  Successive  Line 
Relaxation  scheme. 


C 


APPENDIX  C 


DATA  REDUCTION  PROCEDURE 

In  this  study,  the  bearing  clearance  profile  for  each 
experiment  is  determined  by  fitting  the  experimental  data 
with  an  assumed  surface  profile  using  the  method  of  least 
squares.  The  experimental  data  is  in  the  form  of 
interference  fringes  as  shown  in  Figure  11-13.  Before  the 
fitting  process  can  begin,  the  dark  fringes  on  each  slider 
shoe  from  the  the  two  lasers  are  located  and  their  order 
identified.  To  identify  the  order  of  the  fringes 
appearing  on  the  shoe,  the  location  of  the  fringes  from 
the  Helium-Neon  laser  is  plotted  as  shown  in  Figure  11-14 
by  assuming  their  orders  are  known.  Once  this  is  done, 
the  fringes  from  the  Helium-Cadmium  laser  is  also  plotted 
on  the  same  figure.  If  the  resulting  superposition  of  the 
two  sets  of  fringes  does  not  make  any  sense,  the  process 
is  repeated  until  the  fringes  establish  a  meaningful 
surface  contour  as  shown  by  Figure  11-14.  Once  this  is 
accomplished,  the  order  of  the  "red"  and  "blue"  fringes 
are  determined. 

As  shown  in  Figure  11-14,  in  order  to  determine  the 
flying  profile  of  the  shoe,  that  is,  the  trailing  edge 
clearance  and  the  pitch  angle,  a  "parabolic"  surface 
profile  (crown)  is  fitted  to  the  data  points  on  the  shoe. 
For  example,  in  the  figure,  a  0.051  ym  crown  seems  to  give 
a  good  fit  for  shoe  B  of  the  PA3  slider.  Similar 
procedure  is  followed  for  determining  the  profile  for  shoe 


A  of  the  same  slider.  Because  the  two  shoes  of  the  slider 
are  physically  connected  to  each  other,  their  relative 
positions  with  each  other  must  be  maintained  when  fitting 
the  data.  This  criterion  cannot  not  be  satisfied  when 
fitting  by  hand  and  a  large  error  can  result  in 
determining  the  flying  characteristics  of  the  slider. 

This  potential  error  in  fitting  can  be  significantly 
reduced  with  the  use  of  the  Method  of  Least  Squares 
because  the  flying  profile  determined  using  this  fitting 
technique  is  mathematically  the  best  that  can  be  obtained 
for  the  available  data. 

The  equation  of  the  curve  that  is  to  be  fitted 
through  the  data  points  is  given  by  the  following 
expression:  (refer  to  Figure  11-2  for  notations) 


where 


\  +  hsl  *  V  ■  V 

hf  =  fitted  clearance 


(11.38) 


ht  =  average  trailing  edge  clearance 
0r  =  roll  angle 
Op  =  pitch  angle 

(x,y)  =  location  of  fringe  on  the  shoe 
and  hsl  is  the  surface  profile  given  by  - 


(11.39) 


Therefore,  for  each  fringe  k,  the  location  and  the 
order  of  the  fringe  is  known  and  the  fitted  clearance  is 


SiVL*-i 


c. 


if 


i 


‘■V, 


m 
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given  by  - 


ht  +  hsl„  +  «ryk  "  §pxk 


(11.40) 


Furthermore,  at  the  kjth  fringe,  the  actual  clearance  is 
also  known  and  is  given  by: 

n  *  1,2,3,... 


hfc  ■  (n-1 )— j- 


where 


(11.41) 

n  =  order  of  the  fringe 
X  =  wavelength  of  laser 
hk  s  actual  clearance  at  kth  fringe 

f 

For  each  data  point  in  the  experiment,  there  is  a  h£ 
and  a  h£  associated  with  it.  Because  there  is  error 
between  the  fitted  clearance  and  the  actual  clearance  at 
each  fringe,  it  is  desired  that  the  parameters  ht,  0p,  0r 
be  determined  such  that  the  sum  of  the  squares  of  these 
errors  be  minimized.  Mathematically,  the  problem  can  be 
stated  as  follows: 


Let  the  sum  of  the  squares  of  the  errors  be  given  by 
the  following  expression: 


m 


E  -  i  <hak .  hk7 

k«1 


m  =  no.  of  data  pts.  (11.42) 


and  to  minimize  the  error,  E,  with  respect  to  the 
parameters  ht,  0p,  and  0r  means: 


.  -&L  -  -&L  .  0 

?9r  SPp 


(11.43) 


Substituting  equations  (11.40)  and  (11.41)  into  equation 
(11.42)  and  carrying  out  the  differentiation  as  prescribed 


by  equation  (11.43)*  a  system  of  three  equations  with 
three  unknowns  (ht,  ep,  er)  results: 


The  only  free  parameter  appearing  in  the  above 
equation  is  the  crown  height,  hc,  which  is  in  hsl .  Once 
the  crown  height  is  specified,  the  above  system  of 
equations  (11.44)  can  be  solved  for  ht,  ep,  and  er .  To 
examine  how  well  the  assumed  curve  fitted  the  data,  the 
error  given  by  equation  (11.42)  is  calculated.  Therefore, 
to  determine  the  best  crown  on  each  shoe,  different 
combinations  of  crown  on  each  shoe  are  tried  and  the 
combination  that  gives  the  least  error  in  fitting  is 
designated  as  the  crowns  on  the  two  shoes.  Moreover,  the 

^ti  8  p ,  and  0r  values  calculated  from  equation  (11.44)  for 
the  best  crown  combination  are  also  designated  as  the 
flying  profile  of  the  slider. 
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APPENDIX  D 
I 

ALTERNATING  DIRECTIONAL  IMPLICIT 
AND 

v 

SUCCESSIVE  LINE  RELAXATION  SCHEMES 


This  appendix  gives  an  outline  of  the  Alternating 
Directional  Implicit  (ADI)  scheme  and  the  Alternating 
Directional  Successive  Line  Relaxation  (ADSLR)  scheme  used 
to  solve  the  system  of  equations  resulting  from 
discretizing  the  Reynolds  equation.  To  discretize  the 
equation,  the  finite  difference  approximation  can  be 
applied  directly  to  the  Reynolds  equation  given  by 
equation  (11.25).  It  is  recognized  that  equation  (11.25) 
is  a  mass  balance  equation.  Thus  to  discretize  the 
equation,  a  mass  balance  is  performed  at  each  grid  point. 


C 


C 


Consider  the  variable  grid  mesh  as  shown  below  and 
let  the  arrows  denote  the  directions  of  mass  flow  into  and 
out  of  the  grid  cell  about  node  (ifj)  - 


j-1  j  j+1 


1-1 

L  .  i 

I — 

1 

m—  X 
<□ 

1  V 

i - — 

i  i 

% 

*  X  i  _  , 

i 

1 

i 

1+1 

i 

L ...... . 

1 

-  -  J 

__  AX, 

(*-  AV|- 

AVj 

L  i 

The  mass  balance  at  node  (i,j)  is  simply  stated  as 
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-vT'v  | 

•'.V1'  .  ' 

m, 


follows : 


(mass  in)  -  (mass  out)  =  mass  accumulated 


where  the  mass  fluxes  in  the  x  and  y  directions  are 
given  by  equations  ( 1 1 . 34)  and  (11.35),  respectively. 

Consider  the  mass  in  and  out  of  node  (i,j)  in  the  x 
direction : 


The  mass  flux  into  the  node  is  given  by: 


<*x>in 


C2(Pi-1.j"P»i) 


[e^;(X«-X«-1)  .  l] 


+  C2P i -1 , j 


(11.45) 


while  the  mass  flux  out  of  the  node  is  given  by: 


(mx}out 


C2*PI.l"Pf+1.P 

j"  (r^i+rV 


+  C2Pi ,j 


(11.46) 


Thus  the  mass  balance  in  the  X  direction  can  be 
written  as  follows: 

[A'Pl-1,j  *  A'PIJ  +  C2Pl-l,j  '  A+Plj  +  A+P1+'.j  ‘  C2Pij]  17 


with  'EY  «  i(AYj_^+AYj) 


where 


(11.47) 
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v;%t  ■ 

■ 

V 


A  -  - 


-c2(ap) 


[e^UX)  -  .] 


(11.48) 


with  superscripts  *  snd  *  denoting 
some  averaged  values  between  and  , 

and  Xi  and  Xi  j ,  respectively. 


Similarly  the  mass  balance  in  the  Y  direction  is 
given  as  follows: 

'  B'Pij  +  C4Pl.j-1  *  B+Plj  *  B+PI.J.1  -  <Plj]^ 


with  AX  -  iUXj.j+AXj) 


where 


■c4(aP) 


[.^UY>  -  ,] 


with  superscripts  “  and  +  denoting 
some  averaged  values  between  Yj_1  and  Yj, 
and  Yj  and  Yj+1 .respectively. 


The  mass  accumulated  at  (i,j)  is  given  by: 


(11.49) 


(11.50) 


Hjj  aJCSY 


<1  -  P?.J> 


(11.51) 


Thus,  the  discretized  form  of  the  Reynolds  equation 


becomes : 


[A  P1-1»J  “  A  PU  +  C2P 1-1  ,J  “  A+PIJ  +  A+Pi+1,J  ■  C2Plj] 


■v  '  t 

\ 

It 
i l* 


♦  [B  P1 J-1  "  B  PiJ  *  C4PIJ-1  "  B+Pij  +  B+PI,j+1  "  C4Pij] 


n+1  n 

H,j7»T7(ptJ^r  ptJ) 


(11.52) 


In  this  study,  the  Peaceman-Rachford  ADI  method  is 
employed.  The  advancement  of  the  equation  over  AT  is 
accomplished  in  two  steps  (two  half-time  steps)  - 


1st  Step: 

[n  n  n 

A  PI-1,J  "  A  Pij  +  C2P 1-1 ,  j 

[n+i  _  n+i  _  n+i 

BPI,j-1  “  B  Plj  +  C4PI  ,j«l 


n  n 

A  pij  +  A  Pi+1,j 

.  n+i  .  n+i 
B  PIJ  +  B  Pi,j+1 


n+i  n 
(P,,  -  PfI) 

h!jOT 


2nd  Step: 
n+1 


C2Plj  ]AY 
C4P  5  j  JAX 

(11.53) 


r  .  n+i 

+  [B  PI.J-1 


n+1  _  n+1 

n+1 

n+1 
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It  should  be  noted  that  in  the  first  half-time  step, 
only  the  Y  direction  is  implicit  while  the  X  direction  is 
explicit.  On  the  other  hand,  in  the  second  half-time 


135 


step,  X  direction  becomes  implicit  and  Y  becomes  explicit. 
In  either  step,  the  most  recent  values  are  used  for  the 
explicit  terms;  for  example,  in  the  second  half-time  step, 
the  coefficients  and  the  explicit  terms  are  evaluated  by 
using  the  values  from  the  first  half-time  step. 


I 


C 


c 


c 


c 


c 


It  should  also  be  noted  that  from  examining  equations 
(11.53)  and  (11. 5*0,  it  is  clear  that  terms  involving  P^  j 
from  the  X  and  Y  directions  are  kept  separate  in  the  two 
half-time  steps.  However,  in  Successive  Line  Relaxation 
method,  the  terms  from  the  two  directions  are  assumed 
to  be  implicit.  Therefore,  the  two  half-time  steps  given 


by 

equations 

(11.53) 

and  (11. 
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In  this  study,  the  above  mentioned  scheme  is  called 
the  Alternating  Direction  Successive  Line  Relaxation 
( ADSLR )  scheme.  The  only  difference  between  the  ADI  and 
ADSLR  scheme  is  the  manner  in  which  the  terms  are 

treated . 

Following  the  above  outlined  steps,  two  equations 
similar  to  equations  (11.53)  and  (11.5*0  or  equations 
(11.55)  and  (11.56)  can  be  written  for  each  interior  grid 
point  depending  on  whether  ADI  or  ADSLR  scheme  is  used. 
The  resulting  set  of  equations  when  written  in  matrix  form 
is  tridiagonal  and  can  be  easily  solved  by  using  Gaussian 
Elimination.  Therefore,  to  solve  for  the  steady-state 
solution,  the  scheme  is  iterated  until  the  solution  does 
not  change  between  (full)  time  steps. 


APPENDIX  E 


VARIABLE  GRID  SPACING  ASSIGNMENT  SCHEDULE 

In  this  appendix,  the  scheme  used  to  assign  the 
variable  grid  spacing  in  the  numerical  algorithms  is 
outlined.  Variable  grid  spacing  in  both  the  X  and  Y 
directions  were  used,  but  a  much  finer  mesh  is  used  in  the 
X  direction.  The  same  variable  grid  assignment  scheme  was 
used  for  all  the  different  numerical  algorithms  tested. 
All  the  computer  codes  developed  during  the  course  of  this 
project  can  handle  grid  size  up  to  51  points  along  the 
bearing  and  21  points  across  the  bearing.  However,  for 
all  cases  studied,  21  grid  points  were  used  in  the  Y 
direction  even  when  symmetry  is  present.  Thus  in  effect, 
when  there  is  symmetry,  21  grid  points  are  assigned  across 
half  the  bearing  width  with  the  effect  of  assigning  41 
grid  points  across  the  entire  bearing.  Therefore,  a  more 
accurate  solution  can  be  obtained.  On  the  other  hand,  in 
the  X  direction,  the  number  of  grid  points  assigned  is 
dependent  on  the  operating  bearing  number  based  on  the 
trailing  edge  clearance  with  a  minimum  of  45  grid  points 
and  a  maximum  of  51  grid  points.  Thus,  if  the  bearing 
number  is  very  large,  51  grid  points  would  be  used,  while 
if  the  bearing  number  is  small,  only  45  grid  points  would 
be  used.  The  exact  number  of  grid  points  to  use  is 
decided  by  the  computer  based  on  the  bearing  number.  In 
any  case,  there  would  be  always  9  grid  points  assigned  to 
the  ramp  region.  The  variable  grid  assignment  was 
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designed  such  that  fine  grids  are  joined  with  coarse  grids 
with  grids  of  intermediate  sizes. 

The  scheme  used  to  assigned  the  variable  grid  mesh  is 
given  below  which  is  taken  from  the  computer  program 
developed  during  the  course  of  this  project.  Because  the 
program  was  designed  to  handle  either  constant  grid 
spacing  or  variable  grid  spacing,  the  scheme  used  to 
assign  constant  grid  spacing  is  also  given  below.  It 
should  be  noted  that  unless  the  ramp  length  is  an  integral 
fraction  of  the  "land"  section,  it  would  be  impossible  to 
assign  constant  grid  spacing  throughout  the  length  of  the 
bearing.  Consequently,  there  would  be  one  grid  spacing 
which  will  be  different  from  the  rest.  This  special 
situation  is  fully  accounted  for  in  the  constant  grid 
assignment  scheme. 

The  grid  assignment  scheme  used  for  either  constant 
grid  or  variable  grid  spacing  are  as  follows:  (The 
notation  used  are  self-explanatory) 


c  •••••••AUTOMATIC  GRID  SPACING  ASSIGNMENT 

c  ••*•••••••»•»»•»••••••••»»••»»•»»•••»••*••••' 

c 

C  «»*  GRIDX  -  LOGICAL  VARIABLE 
C  -  .TRUE.  MEANS  X  SPACING  SPECIFIED 

IF(GRIDX)  GO  TO  20 
C  •••  EGRIDX  -  LOGICAL  VARIABLE 
C  -  .TRUE.  MEANS  CONSTANT  GRID 


IF(EGRIDX)  GO  TO  50 

C  •••  VARIABLE  GRID  ASSIGNMENT  IN  X  DIRECTION 
C 

C  NXR  -  NUMBER  OF  DIVISIONS  IN  THE  RAMP  REGION 
C 

NXR=8 

NXR  1=NXR  +  1 
NXL =22 
C 


C  •••  XL  -  LENGTH  OF  THE  BEARING 
C  •••  XLR  -  LENGTH  OF  THE  RAMP  REGION 
C 

DXR=XLR/20. 0 
DXR4=DXR*4 
C 

C  »»•  ASSIGN  GRID  SPACING  IN  THE  RAMP  REGION 
C 

DO  30  1=1,4 
DX(I)=DXR4 
DX(I+4 )sDXR 

30  CONTINUE 
C 

C  •••  ASSIGN  GRID  SPACING  AFTER  THE  BREAK 
C 

DXF =( XL-XLR )/ FLOAT (N XL) 

DXF 14=DXF/4 . 0 
DXFH=DXF/2 . 0 
DXF34=DXF«3. 0/4.0 
DX(NXR 1 )=DXFH 
DX(NXR 1+1 )=DXFH 
DXCNXR 1+2 )=DX34 
DX(NXR1+3)=DX34 
DX(NXR1+4)=DX34 
DX(NXR1+5)=DX34 
NXR6=NXR 1+6 
C 

C  ASSIGN  GRID  SPACING  BEFORE  THE  TRAILING  EDGE 
C 

C  •••  BEARNO  -  BEARING  NUMBER  BASED  ON  THE 
C  •**  TRAILING  EDGE  CLEARANCE 

C  •••  (  =  6*VISC*U*XL/PA/HT/HT) 

C 

TEXL= 1 39.0/15. O/BEARNO 
NXFLAT=NXR  1+NXL+2 
NXFACT=INT(TEXL/DXF) 

NXFLATsNXF  LAT-NXFACT 
TEMP=TEXL-{ FLOAT (NXF  ACT)#DXF) 
ITEMP=INT(TEMP/DXF 14) 

IF0RK=ITEMP+1 

GO  TO  (31,32,33,34),  IFORK 

31  DX(NXFLAT+1 )=DXF 14 
DX(NXFLAT+2 )=DXF 14 
DX(NXFLAT+3 )  =  DXF  14 
DX(NXFLAT+4 )=DXF 14-TEMP 
GO  TO  35 

32  DX(NXFLAT+1 )=DXF 14 
DX(NXFLAT+2 )=DXF 14 
DX(NXF LAT+3)=DXFH-TEMP 
GO  TO  35 

33  DX(NXFLAT+1 )=DXF 14 
DX(NXF LAT+2 )  =  DXF  34-TEMP 
GO  TO  35 


no  ooooo  o  ooo 


34  DX(NXFLAT+1  )sDXF-TEMP 

35  IIsNXFLAT+4-ITEMP+1 
HXFsHXFLAT-6 
DO  36  IsNXR6, NXF 
DX(I )=DXF 

36  CONTINUE 
DX(NXF+1  )  =  DXF34 
DX(NXF+2)  =  DXF34 
DX ( N  XF+3 )  =  DXF  34 
DX(NXF+4)sDXF34 
DX(NXF+5 )=DXFH 
DX(NXF LAT)=DXFH 

•**  ASSIGN  FINE  GRID  SPACING  AT  THE  TRAILING  EDGE 

TEMPsI . O/BEARNO 
DX(II)=3.0«TEMP 
DX(II+1 )=2.0*TEMP 
DX( II+2 )sTEMP 
DX(II+3)=TEMP 
DX(II+4)sTEMP/2. 0 
DX(II+5 )=TEMP/2 . 0 
DX(II+6)=TEMP/4 . 0 
DXCII+7 )= TEMP/4 . 0 
DX ( I 1+8 )= TEMP/5 . 0 
DX ( I 1+9 )= TEMP/5 . 0 
DX(II+1 0)sTEMP/5 • 0 
DX(II+11)sTEMP/10.0 
DX(II+12)=T EM P/15.0 
NX=II+12+1 
NXM1=NX-1 

*•»  END  OF  VARIABLE  GRID  ASSIGNMENT  IN  X  DIRECTION 
GO  TO  20 

•••  CONSTANT  GRID  SPACING  ASSIGNMENT  IN  X  DIRECTION 
•••  NXL  -  NUMBER  OF  DIVISION  WANTED  IN  THE  LAND  SECTION 
50  DDXs(XL-XLR )/NXL 

CHECK  TO  SEE  IF  RAMP  LENGTH  IS  AN  INTEGRAL 
FRACTION  OF  THE  LAND  SECTION 
FACTsXLR/DDX 
NXRsINT(FACT)+1 
C  •••  IF  IT  IS  NOT,  TAKE  CARE  OF  THE  ODD  BALL  WITH  *DX1  • 
DX1 =(FACT-FL0AT(NXR-1 ) )*DDX 
DX( 1 )=DX1 

C  •••  NX  -  TOTAL  NUMBER  OF  GRID  POINTS  (MAX  OF  51 ) 
NXsNXR+NXL+1 
NXR 1 sNXR+1 
NXM1=NX-1 
DO  70  Ia2,NXM1 
DX( I )=DDX 
70  CONTINUE 


oo  o  o  o  o  o  o  o 
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•••  END  OF  CONSTANT  GRID  SPACING  ASSIGNMENT 
•••  IN  X  DIRECTION 

•••  GRID  ASSIGNMENT  IN  Y  DIRECTION 

•••  GRIDY  -  LOGICAL  VARIABLE 

-  .TRUE.  MEANS  Y  GRID  SPECIFIED 
20  IF(GRIDY)  GO  TO  80 
•••  EGRIDY  -  LOGICAL  VARIABLE 

-  .TRUE.  MEANS  CONSTANT  Y  GRID 
IF(EGRIDY)  GO  TO  110 
C  *•*  VARIABLE  GRID  ASSIGNMENT  IN  Y  DIRECTION 
IF( .NOT.SYMM)  GO  TO  81 
C  SYMMETRY  (HALF  BEARING  WIDTH) 

DYSUM=0 . 0 

C  YW  -  BEARING  WIDTH 

DYW4sYW/2. 0/20. 0/M.  0 
DO  91  J=1 ,6 
DY(  J  )  =  DYW4 
DYSUM=DYSUM+DY( J ) 

91  CONTINUE 

DO  92  J=7 , 9 
DY( J )=2. 0*DYW2 
DYSUM=DYSUM+DY( J ) 

92  CONTINUE 

DDYs ( YW/2 . 0-DYSUM) / 11.0 
DO  93  JslO, 20 
DY ( J )=DDY 

93  CONTINUE 

C  END  OF  SYMMETRIC  CASE 

GO  TO  80 

C  NON-SYMMETRIC  CASE  (Y  GRID  SPACING) 

81  CONTINUE 
DYSUM=0 . 0 
DYW4=YW/20. 0/4.0 
DO  94  J  =  1  ,  3 
DY( J )=DYW4 
DYSUM=DYSUM+DY( J ) 

94  CONTINUE 

DO  95  J  =4 , 6 
DY( J)=2. 0*DYW4 
DYSUM=DYSUM+DY( J ) 

95  CONTINUE 

DDY=(Y2/2. O-DYSUM) /4.0 
DO  96  J=7 , 10 
DY( J )sDDY 

96  CONTINUE 

DO  97  Jsl 1 f 20 
J J=2 1 -J 
DY(J)=DY(JJ) 

97  CONTINUE 

C  END  OF  VARIABLE  GRID  SPACING  IN  Y  DIRECTION 
GO  TO  80 
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C 

C  CONSTANT  GRID  SPACING  IN  Y  DIRECTION 
C 

110  CONTINUE 

C  •••  MY  -  NUMBER  OF  GRID  POINTS  IN  Y  (MAX  OF  21) 
MYMIsMY-1 
DDY=YW/MYM1 
DO  120  J=1 ,  MYM 1 
DY( J )=DDY 
120  CONTINUE 

C  END  OF  GRID  ASSIGNMENT  ROUTINE 
80  IF(SYMM)  DY(MY)sDY(MYM1) 

c  mmMmimmmMMiMMiMmmmiiiii 
C  mimmMmmMmmMiMimimmmi* 
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